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crS i Optical and near-infrared variability is a well-known property of young stellar 

Tij- ! objects. However, a growing number of recent studies claim that a considerable 

^ ! fraction of them also exhibit mid-infrared flux changes. With the aim of studying 
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and interpreting variability on a decadal timescale, here we present a mid- infrared 
spectral atlas containing observations of 68 low- and intermediate mass young 
stellar objects. The atlas consists of 2.5-11.6 /im. low-resolution spectra obtained 
with the ISOPHOT-S instrument on-board the Infrared Space Observatory (ISO) 
between 1996 and 1998, as well as 5.2-14.5 /im low-resolution spectra obtained 
with the IRS instrument on-board the Spitzer Space Telescope between 2004 and 
2007. The observations were retrieved from the ISO and Spitzer archives and 
were post-processed interactively by our own routines. For those 47 objects where 
multi-epoch spectra were available, we analyze mid-infrared spectral variability 
on annual and/or decadal timescales. We identify 37 variable candidate sources. 
Many stars show wavelength-independent flux changes, possibly due to variable 
accretion rate. In several systems, all exhibiting 10 fim silicate emission, the 
variability of the 6— 8/xm continuum and the silicate feature exhibit different 
amplitudes. A possible explanation is variable shadowing of the silicate emitting 
region by an inner disk structure of changing height or extra silicate emission 
from dust clouds in the disk atmosphere. Our results suggest that mid-infrared 
variability, in particular the wavelength-dependent changes, are more ubiquitous 
than was known before. Interpreting this variability is a new possibility to explore 
the structure of the disk and its dynamical processes. 

Subject headings: stars: pre-main sequence — stars: circumstellar matter — 
infrared: stars — methods: data analysis — techniques: spectroscopic 



1. Introduction 

At the early phases of their formation, stars arc intimately linked to their circumstellar 
environment. Dust particles in the circumstellar disk and envelope are heated both internally 
by viscous energy released in the midplane of the accretion disk, and externally by irradiation 
from the central star. The absorbed energy is re-radiated as thermal emission at infrared (IR) 
wavelengths. The shape of the IR spectral energy distribution (SED) reflects the temperature 
and density structure of the circumstellar environment. Of particular importance is the 



*This work is based on observations made with the Infrared Space Observatory (ISO) and with the Spitzer 
Space Telescope. ISO is an ESA project with instruments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the UK) and with the participation of ISAS and NASA. 
Spitzer is operated by the Jet Propulsion Laboratory, California Institute Technology under a contract with 
NASA. 
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mid-IR domain (3-25 /xm), where the disk emission originates from the region where planets 
are formed. The same domain contains spectral features of PAH particles, small silicate 
grains, and molecular ices. PAH features gi ve infornaation about the UV radiation field and 
ionization balance in the circumstellar gas ( iTielensI 120081 ) : the shape and amplitude of the 
10 /im silicate feature is partly related to the degree of crystallinity of the silicate grains, 
indicating how far these grains evolved from the interstellar medium (jHenningj 120101 ): while 
ice featur es are used to study the organic chemistry and thermal history of the circumstellar 
material (Ivan Dishoeckll2004j ). 



Although optical and near-IR variability is a well-known property of young stellar ob- 
jects (YSOs), a growing number of recent stud ies claim that a cons i derable fraction of them 



also exhibi t mid-IR photometric changes (e.g. iBarsony et al.l 120051 : [Morales- Calderon et al. 



20091 . 12OIII ). In young eruptive stars and some TTauri stars, changing accretion is thought to 

be re sponsible for mid-IR variability (e.g. [Abraham et al. iKospal et al. Ivan Boekel et al, 
20101). There are claims that in some young eruptive stars the temporarily increased heat 
leads to structural changes in the innermost part of the system, resulting in decreased ex- 



tinctio n and consequ ently in IR brigh tness changes (see, e.g., V16470ri in iMuzeroUe et al. 



2005bl . or PVCep in iKun et al.ll201ll ). In some intermediate mass young sta rs, structural 



changes in the iiiner d isk were invoked to explain their mid-IR flux changes (jJuhasz et al. 



20071 : ISitko et al.ll2008[ ). A similar reason, i.e. changing inner disk height, r aay also explain 



the yariability of those low -mass objects whose disks contain holes and gaps (IMuzeroUe et al. 



2OO9I : lEspaillat et al.l 1201 if ). These indications of dynamical processes suggest that disks are 
less steady-state than thought before. 

The wavelength-dependence of the mid-IR variability carries information on the disk 
and envelope structure, because the emission at different wavelengths corresponds to disk 
regions with different temperatures and thus to different radial distances from the star. An 
interesting aspect is the exploration of possible changes in the strength and profile of IR 
spectral features. Several studies published so far concentrate on the variability of the 10 /im 
silicate feature. Time-dependent shadowing of disk areas by the in ner gas disk could explain 
the spectral changes in the Herbig stars HD 31648 an d HP 163296 (jSitko et al.ll2008l) . In the 
highly accreting TTauri stars DGTau and XZTau ( iBary et al.l 120091 ). self-absorption may 
play a role in the variability of the 10 ^m feature. Changes in the minera logical composition 
of the emitting dust particles were also observed by lAbraham et al.l ( 120091 ) in the most recent 
outburst of EX Lupi, where they detected the appearance of crystalline silicates that were 
not present in quiescence. However, more detailed investigations are limited by the lack of 
multi-epoch mid-IR spectroscopic databases. 



The timescale of mid-IR variability is related to the timescale on which the heating of 
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the emitting region can vary and on which the dust particles in the disk can react to changes 
in the energy input. For example, almost instantaneous (seconds to minutes) responses to 
fluctuations in the stellar illumination may indicate emission from an o ptically thin medium, 
while longer timescales (years) may indicate an optically thick medium ( IChiang &: Goldreich 
19971 ). Dynamical changes leading to rearrangement of the density structure can occur on the 
timescale of an orbital period. Thus, the timescale of the variability is an essential input to 
understand the physics behind the flux changes. Recent studies carried out with the Spitzer 
Space Telescope sampled MIR variabilit y on several different tirnescale s. While hourly vari- 
ations are rare and small in amplitude (IMorales-Calderon et al.ll2009l ). variability on daily, 
weekly, and yearl y timescales is ubiquitous, and amplitudes rarige from a few tenths of magni- 



tudes to >lmag ([Morales- Calderon et al.ll2009l : iLuhman et al 



2010: Morales-Calderon et al. 



201 ll : lRebullll201l[ ). These timescales, if caused by changes in the disk structure, indicate 
that the variability is connected to the innermost 0.1-1 AU of the disk. With Spitzer alone, 
however, one cannot investigate whether variability on longer timescale exists. Such data 
would provide information on dynamical changes at larger radial distances from the star, 
typically several AU. Such study can only be performed by comparing data from two IR 
space missions separated in time by several years. 

Motivated by the above mentioned considerations, here we present a new atlas of 2.5- 
11.6 /im low-resolution spectra obtained with the ISOPHOT-S instrument on-board the In- 
frared Space Observatory (ISO) between 1996 and 1998, as well as 5.2-14.5 fim low-resolution 
spectra obtained with the IRS instrument on-board the Spitzer Space Telescope between 2004 
and 2007. The overlapping wavelength range and different epochs of ISO and Spitzer make it 
possible to do a systematic mid-IR spectral variability study for a large sample of well-known 
YSOs both on annual timescale (less than a few years, i.e. within the cryogenic lifetime of 
either telescope), and of nearly a decade (by comparing the results of the two telescopes). 
Although the beams of the two instruments signiflcantly differ, beam confusion is usually 
not an issue, because at mid-IR wavelengths, the size of the emitting region (warm dust in 
the inner few Alls) is smaller than either in the optical (scattered light) or in the far-IR/mm 
(thermal emission from cold dust). We identify typical variability patterns and hypothesize 
about their possible physical origin. Our results on monitoring and interpreting variability 
may provide a powerful "extra dimension" of information on the structure of the circum- 
stellar material. Our sample may also serve as a starting point to plan follow-up mid-IR 
monitoring campaigns. 
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2. Observations and data reduction 

2.1. Sample selection 

We made a query of the ISO data archive for ISOPHOT-S observations (AOT PHT40) 
and selected observations of low- and intermediate mass YSOs by consulting the SIMBAD 
database and the literature. The query resulted in 94 measurements of 68 YSOs; including 
35 low-mass and 33 intermediate mass stars. This forms the target hst of the present study. 
Then we checked the Spitzer archive for 5.2-14.5 fj,m low-resolution IRS spectra. We found 
that such Spitzer spectra were available for 51 of our targets. Table [1] presents the list of our 
targets and the log of observations. 



2.2. ISO data reduction 

ISOPHOT-S was a sub-instrument of ISOPHOT on-board the ISO. It produced low- 
resolution mid-IR spectra of 128 data points covering the 2.5-4.9 (ISOPHOT-SS) and 5.8- 
11.6 nm (ISOPHOT-SL) wavelength ranges with a spectral resolution of R ^100 during the 
active period of ISO (1995-1998). The aperture was 24"x24". For a detailed description 



of the ISOPHOT-S instrument and its calibration strategies, see lAcosta-Pulido fc Abraham 



(]2003f ). 



Most observations presented here consist of pairs of separate ON/OFF measurements, 
but chopping between the source and background position(s), as well as small scans/maps 



were also utilized. Edited Raw Data 



software (OLP v. lO.O. lLaureijs et al. 



ERD ), produced by the ISOPHOT Off-Line Processing 



20031 ). were downloaded from the ISO Data Archive for 



all of our targets. Standard processing to Auto Analysis Product Data (AAP) level was per- 



forme d with the ISOPHOT Interactive Analysis Software Package (PIA vlO.O. iGabriel et al. 



19971 ). Further processing was done using self- developed IDL routines consisting of the fol- 
lowing steps: improved deglitching, correction for memory effects due to preceding bright 
source measurements, background subtraction in the lack of dedicated OFF observation, 
correction for flux loss in off-center pointings (in case the offset was larger than 1"), and 
an empirical photometric correction. Several of these algorithms were aimed at correcting 
artifacts that were discovered after the development of the ISO offline processing software 
was finished. We found that in a number of cases, the reprocessing significantly improved 
the data quality. According to our systematic study of 28 normal star spectra, we concluded 
that the uncertainty of the absolute calibration is 5-10% for bright sources (brighter than 
a few Jy) and 30— lOOmJy for fainter sources (<1 Jy). Note that these figures are valid for 
each independently calibrated pixel. A detailed explanation of the correction steps and the 
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error budget can be found in Appendix A. 



2.3. Spitzer data reduction 

The Infrared Spectrograph (IRS) on-board the Spitzer Space Telescope had two low 
resolution and two high resolution modules, covering the 5.2-38 //m wavelength range. Since 
our aim is to look for mid-IR variability, here we utilize the Short-Low channel that covers 
the 5.2-14.5 //m wavelength range, thus has the largest wavelength overlap with ISOPHOT- 
S. The Short-Low channel has a spectral resolution of R ~ 60-127, comparable to that of 
ISOPHOT-S. The width and length of the IRS Short-Low slit is 3'.'6 and 57", respectively. 
According to the IRS instrument handbook (Version 4.0), the uncertainty of the absolute 
flux calibration is 5-10%, while the repeatability is 2-5%. 

We downloaded basic calibrated data (BCD) files from the Spitzer archive, processed 
with the pipeline version 18.18.0. Measurements were taken either in staring mode (two 
nod positions), or in spectral mapping mode (small 2x3 maps). In the latter case, we 
used the central 2 positions, and considered them as normal staring measurements. We 
first subtracted the two nod positions from each other. Then we used the Spitzer IRS 
Custom Extraction software (SPICE) to extract the positive signal from the two-dimensional 
dispersed images. We extracted the spectra from a wavelength-dependent, tapered aperture 
around the target, and we averaged the two spectra corresponding to the two nod positions. 
In case the target was not well-centered in the slit (i.e. offset from the slit center perpendicular 
to the slit was larger than 0.5"), we corrected the sp ectra for flux loss using the measured 



IRS beam profiles (for details, see lKospal et al.ll2008l ). 



2.4. Comparison of the ISOPHOT-S and Spitzer/IRS calibration 



The absolute spectrophotometric calibration o 
cal spectral templat es of standard stars derived by 



ISOPHOT-S was based both on empiri- 



Cohen (120031) . and on stellar atmosphere 



models provided by iHammersley fc Jourdain de Muizoru ( l2003l ) using infrared observations. 
Spitzer/IRS spectra, on the other hand, were c alibrated using standard stars for which model 
atmospheres are available in lDecin et al.n 20041) . Nevertheless, Cohen templates were used to 
verify the calibration of IRS (iHouck et al.ll2004j ). Thus, despite the differences in calibration 
strategy, we do not expect significant systematic discrepancies between the ISOPHOT-S and 
the IRS spectra. In order to check this, we reduced the spectra of HD 181597 (HR7341), a 
Kl HI type star that was used as a primary low resolution standard for Spitzer and was also 
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used in the calibration of ISOPHOT-S. In Fig. [T]we plotted the observations along with the 
template and model spectra, and we also calculated their ratio. The figure shows that the 
template and model spectra are close to each other, and the observations are within 10% (for 
ISOPHOT-S) or 5% (for Spitzer/IRS) of the template/model spectra. We emphasize that 
no wavelength-dependent trend can be seen in the ratios, suggesting that the calibration 
uncertainties would not affect the observed colors of the targets. Thus, we conclude that the 
accuracies of the absolute calibration mentioned in Sec. 12.21 and 12.31 are plausible and there 
is no need to increase the error bars when comparing spectra obtained with the two different 
telescopes. 



3. Catalog of observations and the spectral atlas 



In Table [T] we present the catalog of spectroscopic observations. The table contains the 
following columns: 

Column (1): Name of the target as common in the literature. 

Column (2): L - low-mass pre-main sequence stars (T Tauri-type stars and embedded young 
stellar objects); I - intermediate mass pre-main sequence stars (Herbig Ae/Be stars). B, A, 
and F stars are considered intermediate mass stars, while G, K, and M stars are considered 
low-mass stars. Spectral types or stellar masses for sou rces where SIMBAD d oes not pro 
vide s pectral types were t aken from IStraizvs et al.l ( 120021 ) for IR AS 03260-1-311 l.lMyers et al. 



(119871 ) for BAR N 5 IRS 1. iBrinch et al. 



LDN 15 51 IRS 5, 



Sandel 



WL 16 , lBontemps et al 



fc Weintraub 



(2007) for LDN 1489 IRS, iLim fc Takakuwa 



mm forReiDurth50 N IRS 1 



Simon et al. 



mm for WL 6, iHodapp et al.l (Il996[ l for 00 Ser, and 



( 2006 ) for 
(1995) for 



Pontoppidan et al 



mm for [SVS76]Ser4. 



Column (3): Spectral types from the SIMBAD database. 

Column U)-(5): J2000 co ordinates of the source from the 2MASS All-Sky Catalog of Point 
Sources ( ICutri et al.ll2003l ). except where indicated in the table. Since the peak of the optical 
and IR emission might not coincide, we decided to use coordinates from 2MASS instead of 
optical positions. These coordinates are precise enough to carry out the offset-correction if 
the source was off-center in the ISOPHOT-S aperture or in the Spitzer/IRS slit. 

Column (6): Unique 8-digit identifier of the on-source ISOPHOT-S observation in the ISO 
Data Archive (TDT number). 

Column (7): Corrections applied to the ISOPHOT-S observations: M - memory correction; 
B - in the lack of dedicated OFF measurement, zodiacal background was predicted and 
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subtracted; O - position offset correction applied. An exclamation mark signals spectra 
where the uncertainty of the correction is above 10%, thus they are not used for variability 
analysis (for details, see Appendices A and B). 

Column (8): Date of the ISOPHOT-S observation. 

Column (9): Unique identifier of the Spitzer/IRS observation in the Spitzer Data Archive 
(AOR). 

Column (10): Corrections applied to the Spitzer/IRS observations: O - position offset 
correction applied. 

Column (11): Date of the Spitzer/IRS observation. 

Column ( 12)-(14 ): Inventory of clearly visible spectral features. Si-0 - 10 /im silicate feature; 
Ices - H2O 3.1 ^m, CO2 4.27 //m, H2O 6.0 /im, CH3OH 6.8 yum, or CH4 7.7 /xm; PAH - 3.3, 
6.2, 7.7, 8.6, 11.2, or 12.7 /xm bands of polycyclic aromatic hydrocarbons, em - feature in 
emission; abs - feature in absorption; . . . - feature not present. 

Column (15): Type of the spectral shape: PAH dom. - the spectrum is dominated by 
PAH emission features, no 10 //m silicate feature is discernible; sil. em. - the 10 /xm silicate 
feature is in emission; sil. abs. - the 10 /xm silicate or ice bands are visible in absorption; 
. . . - undecided; (for details, see Sec. H]). 

Column (16): Variability: yes - candidate variable source; no - constant source; . . . - unde- 
cided. 

The table and the spectra are also available electronically. 

4. Results 

In Fig. |2]we present the ISOPHOT-S and Spitzer/IRS spectra for each object. Where 
clearly detected, we marked with vertical dashed lines the wavelengths of molecular ice 
bands (H2O 3.1 /im, CO2 4.27 /xm, H2O 6.0 /xm, CH3OH 6.8 /xm, CH4 7.7 /xm); and with 
vertical dotted lines the wavelengths of PAH features (3.3, 6.2, 7.7, 8.6, 11.2, 12.7/xm). The 
figures show both cases where the plotted spectra agree, and sources where the measurements 
at different epochs seem to differ. Some of the differences can be explained by serious 
instrumental artifacts or source confusion resulting from the different beam sizes of the 
ISOPHOT-S and IRS instruments. We studied the literature of each object in order to 
establish whether it is part of a multiple stellar system and whether there is any indication 
for extended near or mid-IR emission in the vicinity. These information are presented in 
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Appendix B for the affected objects. Based on tliis, we could decide wlietfier it is meaningful 
to compare the different spectra in Fig. [2]to look for temporal variability. For certain objects, 
Appendix B also contains a non-exhaustive comparison with other mid-IR spectra published 
in the literature. 

In order to explore the significance of the observed flux changes and check their wavelength- 
dependence, we computed the ratio of spectra taken at different epochs. In Fig. [3]we plotted 
these flux ratios for sources where spectra at more than one epoch were available and their 
comparison was meaningful (47 stars). Where several spectra existed, we took the two most 
extremes. In case the ratio of a Spitzer/IRS and an ISOPHOT-S spectrum is plotted, the 
former was first resampled to the latter's wavelength resolution. Error bars plotted in this 
figure represent the quadratic sum of the errors of each spectrum. Where the ratio of two 
Spitzer spectra is computed, instead of the 10% absolute calibration uncertainty, we used 
the 5% uncertainty of the repeatability (BFOri and HD 98800). We also overplotted with 
vertical dotted/dashed lines the wavelengths of PAH/ice features for sources where these 
features are present. 

For the purpose of further analyses, we classified our sources into three groups based on 
the shape of the mid-IR SED and the presence of various spectral features as follows: 

• PAH dominated sources: the spectrum is dominated by PAH bands, no apparent 10 /im 
silicate emission band is visible (13 objects); 

• Silicate emission sources: the spectrum is dominated by the 10 /xm silicate emission 
band, otherwise it is rather featureless (36 objects). In a few sources, weak PAH 
emission is present; 

• Silicate absorption sources: the spectrum is rising towards longer wavelengths, it dis- 
plays 10 /um silicate absorption and/or several ice absorption bands (16 objects). 

Our classification is presented in Column (15) of Tabled! as well as in Fig. |2]in parentheses 
under the name of each source. We could classify all sources except DGTau, Cedll2IRS4, 
and CKl, which display an apparent combination of silicate emission and absorption. 

4.1. Statistics and wavelength-dependence of mid-IR variability 

Inspecting Figs. |2] and |3l we could conclude that the flux ratio in many sources differ 
from unity at a level of more than la, for part of or the whole wavelength range. Note that 
la here refers to individual wavelength channels, and averaging the spectrum for a wider 
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wavelength range suppresses the random, channel-to-channel noise. Thus, a 10-15% flux 
change for a continuous wavelength range is in fact already a strong indication of variability. 
Out of the full sample of 68 stars, variability could not be decided for 21 objects, for having 
only one spectrum, or due to problems related to faintness, source confusion, or instrumental 
artifacts. We found 10 stars (21%) to be constant between 6 and 11.5 /im. The remaining 37 
sources (78%) were classified as variable objects. These are marked in Column 16 of Tabled] 

The wavelength-dependence of the observed variability can also be characterized from 
our data. By looking at Figs. [2] and El we conclude that the wavelength-dependence of the 
flux changes is mostly smooth over the 5-12 fim spectral range, and radical changes within a 
few channel distance are not observed. Although changes in the strength of the 10 ^m silicate 
emission feature are ubiquitous and will be quantified and discussed later (see Sections 14.31 
and 15. 2p . we found no sources where the silicate or any other spectral feature emerged or 
completely disappeared from one epoch to the other. Neither did we find sources where the 
shape of the feature changed significantly, implying that alteration of the size distribution 
or mineralogical cor nposition of the emitt ing dust grains (cf. the crystallization during the 
outburst of EXLup, [Abraham et al.ll2009l ) is not widespread in our sample. 



4.2. Synthetic photometry and variability timescales 

The smooth trends of the flux ratio curves in Fig. [3] suggest that the wavelength- 
dependence of the variability can be fully described, with high signal-to-noise ratio, via 
integrating the spectra over a few wavelength intervals and looking at the changes of these 
integrals. Following this idea, we derived synthetic photometry from the spectra. For objects 
with silicate emission or absorption, we calculated synthetic photometry for each source by 
averaging the ratio curves in Fig. |3] between 6 and 8 ^m and between 8 and 11.5 /xm, and 
converting them to magnitude differences. The former wavelength range is representative of 
the continuum, while the latter one covers the 10 fim silicate feature and is close to the usual 
N photometric band. For PAH dominated objects, we defined a "continuum band" by aver- 
aging ratios at 6.5-7 and 9.5-10.5 fim and a "PAH feature band" containing the remaining 
parts of the spectrum. 

The synthetic magnitudes can be used to assess quantitatively a fundamental ques- 
tion: whether the characteristics (frequency and amplitude) of the variability are the same 
on annual timescales (ISOPHOT-S vs. ISOPHOT-S or Spitzer vs. Spitzer) and on decadal 
timescales (ISOPHOT-S vs. Spitzer), or there is a hint for long-term variability trends (higher 
frequency or larger amplitude) not measurable by ISOPHOT-S alone or Spitzer alone. In 
Fig. m we compared the distribution of magnitude changes in the annual and decadal cases 
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for the three objects types defined above. Histograms are plotted separately for the "con- 
tinuum" and the "feature" magnitude changes. For PAH dominated objects, no systematic 
difference can be seen between the short-term and the long-term distributions, either in 
the continuum or the feature (Fig. H^, e). For objects with silicate emission, we plotted 
separately the intermediate mass and the low-mass stars, but even with this separation, no 
significant differences can be seen between the two timescales (Fig. |Dd, c, f, and g). The only 
weak hint for such long-term trends might be seen for silicate absorption sources, where both 
in the continuum and the feature, the decadal changes seem to be somewhat more common 
and exhibit somewhat larger magnitude differences than on annual scale (Fig. HJi, h). Al- 
though the number statistics are limited, we can conclude that the characteristic timescales 
of the mid-IR variability we are observing is less than or equal to a few years, and in the 
following we do not differentiate between annual and decadal changes when studying the 
physical reasons of the variability. 



4.3. Trends in the mid-IR flux changes 

In Fig. [5] we plot the obtained synthetic magnitude changes. As in Fig. [3], we always 
considered the two most extreme spectra. Sources with constant fiux in the whole 6-11.5 fim 
wavelength range (Fig. E]) are situated close to the (0,0) point. Most of them fall into an 
area within 0.1 mag, indicated by a dark gray circle. Moreover, we checked the synthetic 
photometry for the calibration star HD 181597, which was observed by both ISOPHOT-S and 
Spitzer/IRS (Fig. [1]), and found Amag(6-8 /im)=0.067mag, Amag(8-11.5 /im)=0.076 mag. 
For these reasons, we consider the value of 0.1 mag as a practical threshold between variable 
and constant sources. Objects with a wavelength independent fiux change (either brightening 
or fading) are distributed along the 45 degree diagonal line (for guiding the eye, this line, 
broadened by a ±0.1 mag typical uncertainty, is plotted with light gray). A similar, but 
vertical stripe is plotted for objects with silicate emission, marking the location of those 
sources where the continuum part was constant, but the 10 silicate emission feature 
varied significantly. 

Several trends in the distribution of the points can be recognized in Fig. [51 

• Objects whose spectra dominated by PAH emission are plotted in the upper left panel 
of Fig. O Out of these 8 objects, 5 are marked as constant in column 15 of Tab. [TJ 
These are situated within 0.1 mag of the (0,0) point. All the remaining three sources 
(RRTau, CUCha, and HD 135344 B) are distributed along the diagonal stripe with 
a magnitude difference of 0.2-0.3 mag both in the continuum and in the PAH feature 
band. We suggest that these objects show real physical variability (RRTau exhibited 
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flux changes of 0.5 mag also between 2.5-5 /im, see Fig. [2]). Tlie variability in the 6- 
11.5 yum range seems to be wavelength independent, as indicated by the distribution of 
points along the diagonal stripe. This means that for each source the variation in the 
PAH features and in the continuum has similar amplitude. 

• Stars with silicate emission are plotted in the upper right (Herbig stars) and lower 
left (T Tauri stars) panels of Fig. [51 They exhibit the largest diversity in variability 
patterns within the whole sample. Three Herbig stars (HD 95881, HD 142527, and 
AKSco) and no T Tauri stars are constant. 11 stars show wavelength independent 
flux changes of >0.1mag, located along the diagonal stripe (Herbig stars: BFOri, 
HD 104237, HD 144432, VVSer, WWVul; T Tauri stars: DRTau, VZ Cha, WXCha, 
WWCha, EXLup, SCrA). The amplitude can be as high as 0.7 mag. Remarkably, 
another 11 sources change their brightness primarily in the 8-11.5 /im band (amplitude 
is <0.6mag) while at shorter wavelengths, they are less variable or constant (Herbig 
stars: HD 98800, HD 139614, HD 142666, 51 Oph; T Tauri stars: UZTau, VYTau, 
CRCha, CTCha, VWCha, Glass I, CVCha). Interestingly, one T Tauri star does 
not follow any of the above described trends. XX Cha became brighter at shorter 
wavelengths and fainter at longer wavelengths, exhibiting an anti-correlation between 
the two synthetic photometric bands. Although the amplitude of the variability of this 
source is small, the phenomenon is probably real, since it is the slope of the mid-IR 
SED that changed, rather than the more uncertain absolute flux level only. 

• Stars with silicate and/or ice absorption (lower right panel in Fig. [5]) show frequent vari- 
ability (only one source out of 10, Ced 111 IRS 5, is constant). The observed variations 
mostly follow the diagonal stripe, indicating wavelength independent flux changes. The 
amplitude of the variability is at most 0.4 mag. Slight deviations from this trend can 
be seen in those objects where the 10 fim silicate absorption feature is very deep, thus, 
their 8-11.5 //m flux is close to zero and remains essentially constant (e.g. Haro6-10, 
LDN1551IRS5, OOSer). 

There are objects where the detected flux changes are only marginal (>0.1mag), and 
their association with the above mentioned trends may not be unambiguous. However, for 
each trend, there are clear examples with magnitude changes >0.3mag, which proves the 
existence of these trends. The association of a certain source with a variability trend might 
depend on which measurements are compared. E.g. in the case of WW Cha, we have one 
ISOPHOT-S and three Spitzer/IRS spectra. The comparison of the Spitzer spectra indi- 
cate wavelength-independent changes, while the ratio of the ISOPHOT-S spectrum and the 
Spitzer/IRS spectrum obtained on 07-Mar-2006 provides a clear example where the flux 
variations above and below ^9 fim anti-correlate, similarly to XX Cha. VYTau, BFOri, 
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CR Cha, and Glass I are further examples for cases where certain spectral ratios show wave- 
length independent flux changes while others indicate changes only in the silicate feature 
or hints for anti-correlation between the continuum and the feature. This suggests that the 
trend of variability may not be a unique characteristic of a given source, but may change in 
time. 

For 15 sources where two or more ISOPHOT-S spectra were available, we could also 
study variability in the 2.5-4.9 /im wavelength range. We found that five sources (UZTau, 
VYTau, HD 97300, CVCha, HD 135344 B) are constant between 2.5-4.9 /xm. Out of these, 
UZTau, VYTau, and CVCha are slightly variable at longer wavelengths, while HD 98300 
and HD 135344 B is constant throughout the whole 2.5—11.6 fim domain. The two variable 
PAH dominated objects in this sample (RRTau and CUCha) have larger flux changes at 
shorter wavelengths than at A > 5 //m (0.2-0.5 mag). A number of silicate emission objects 
show variability amplitudes comparable to those at longer wavelengths (DRTau, Glass I, 
HD 104237, and EXLup). One silicate emission object, however, displays anti-correlation 
between the shorter and longer wavelength regimes (CRCha). This behavior is similar to 
what we observed for WW Cha and XX Cha at longer wavelength, except that the constant 
pivot point is at somewhat shorter wavelength. The three variable silicate absorption objects 
(Haro6-10, HH 100 IRS, and 00 Ser) display slightly higher variability amplitude in the 2- 
5 /im range than at longer wavelengths. 



5. Discussion 



In the previous sections we investigated the characteristics of variability for a sample 
of well-known YSOs. Based on the magnitude changes plotted in Figs. H] and out of 
47 objects where multi-epoch observations were available, we found 37 stars (79%) to vary 
> 0.1 mag and 20 stars (43%) to be variable at > 0.3 mag level either in the continuum or the 
feature. Our results can be compared with those of ear lier mid-infrared vari ability studies 
carried out on annual (typically 1-4 years) timescales. iBarsony et al. surveyed the 

pOph cloud and found that at least ~20% of the YSOs are v ariable at 10 /xm on a few years 
timescale. Based on multi-epoch Spitzer/IRAC photometry, iLuhman et al.l (120101 ) reported 
that 44% of Class O/I/II sources in the Taurus, and 26% in the ChamaeleonI star forming 
regions exhibit at le ast 0.05 mag flux chan ges. Recently, from the analysis of multi-epoch 
Spitzer/IRS spectra, lEspaillat et al.l ( 120111 ) concluded that 86% of the transitional and pre- 
transitional disks in the Taurus and Chamaeleon star forming regions display larger than 
10% flux variability in the 5-38 /xm wavelength range. Our results confirm the high incidence 
rates found in these studies and demonstrate that in our - rather heterogeneous - sample of 
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young objects, annual/decadal mid-IR variability is ubiquitous. 



Our data set also offers a possibility to explore the characteristic timescales of the 
processes responsible for the observed variability. These timescales could place constraints 
on the nature of the underlying physics. In Section 14.21 and Fig. H] we statistically proved 
that the distribution of the variability amplitudes on annual and on decadal timescales 
do not differ significantly. This new result implies that on decadal timescales we observe 
the same variability process than on shorter timescales, and no extra long-term physical 
mechanism is evident from our data. This fact constrains the physical timescale of the 
variability to less than a few years. There are hints that this characteristic timescale may 
even be shorter. In our Spitzer/IRAC monitoring program (PID 60167, PI: P. Abraham), 
we obtained 14-day-long light curves at 3.6 and 4.5 /xm of 38 selected low- and intermediate 
mass YSOs and found that practicall y all targets were variabl e with amplitudes from a few 
tens of millimagnitudes to 0.4 mag. [Morales- Calderon et al.l (120091 ) observed 69 YSOs in 
the IC 1396A dark globule between 3.6-8.0 with Spitzer/IRAC. The obtained 14-day- 
long light curves indicated that more than half of the observed YSOs show daily variations, 
with amplitudes ra nging from 0.05 mag to 0.2 mag. Finally, th e initial results of the Spitzer 
YSOVAR program (lRebullll201ll : iMorales-Calderon et al.ll201ll ). based on the analysis of 40- 
day-long light curves obtained for sources in the Orion Nebula Cluster and other well-known 
star-forming regions, indicate that 70% of the YSOs with IR excess are variable with days to 
weeks timescales and with amplitudes up to higher than 1 mag. These papers suggest that 
the incidence and the amplitude of flux changes we found on decadal timescale seem to be 
present already on the much shorter timescale of a few weeks. If confirmed, this result would 
strongly constrain the characteristic timescale of the physical variability process, possibly 
pointing to the innermost part of the system as the region of the true cause of variability. 

In order to examine flux changes, we integrated our spectra in relatively broad synthetic 
photometric bands. Consequently, the amplitude of variability within limited wavelength 
ranges can be even higher. Moreover, since our conclusions are based on the comparison of 
only two or three different spectra, the incidence numbers are lower limits, and the fraction 
of mid-IR variable stars may very well be larger. A dedicated mid-IR monitoring program 
would probably reveal variability for many of those objects that we observed to be constant. 
The amplitude of variability might also turn out to be larger if more epochs are taken 
into account. Examples are HD 169142 in Fig. ^ HD 104237 in Fig. M and HP 1 42527 
in F ig. [T71 where we plotte d additional mid-infrared spectra from ISylvester et al.l (119961 ) 
and Ivan Boekel et al.l (120051 ). Based on these results one might predict that a very high 
fraction of all low- and intermediate mass young stars are variable at a certain level in the 
mid-IR wavelength regime. Nevertheless, our data also revealed significant differences in the 
incidence of variability among objects with different SED types. As an example. Figures H] 
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and [5] clearly show that the fraction of variable stars is lower among the PAH-dominated 
objects than among the TTauri stars. In the following we discuss the possible physical 
reason of variability for each object type separately. 



5.1. PAH dominated objects 



In our sample, sources whose spectra are dominated by PAH emission are all interme- 
diate mass stars with spectral types ranging from BO to F6, except the GO-type TTauri 
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from molecules transiently heated by the UV photons o f the central source. Thes e species 
can be located at relatively large distances from the star. iHabart et al.l (12004 l2006l ) resolved 
the PAH emission around several Herbig stars (e.g. HD 169142, CUCha, HD 100453, and 
HD 100546) and found that PAH emission is originated from the illuminated surfaces of flared 
disks. The 3.3 /im feature is the least extended {R < 30 AU), while the 6.2, 7.7, and 11.3 /im 
features are emitted by the outer disk region (i?~100 AU). In other cases, the source of 



PAH emission is an even more extended envelope (e.g. HD 97300. ISiebenmorgen et al.lll998 



Kospal et al.ll2012l ). or mig ht be coming from the surrounding cl oud material (e.g. 
03260-1-3111 and MWC865. Isioan et al.l[l999l : ISoersma et al.ll2009[ ). 



IRAS 



When explaining the spectral variability of PAH dominated objects, one has to take into 
account that the spatial location of the regions emitting the mid-IR continuum may differ 
from that of the PAH emission. The bulk of the continuum emission at mid-IR wavelengths is 
radiated by warm, large dust grains in thermal equilibrium. Depending on the temperature 
and luminosity of the central star, these grains are conflned to the inner disk, within a 
few Alls from the star, much closer than the source of the PAH emission. Our spectra 
do not exhibit any discernible silicate emission around 10 fim. In a few cases, the PAH 
emission is so strong that a weak silicate emission might be present but difficult to notice 
(IRAS 03260+3111, CUCha, HD 97300, WL 16, MWC865). In other cases, the underlying 
continuum can be clearly seen and the silicate emission is apparently absent (HD 34700, 
RR Tau, HD 100453, H D 135344 B, HD 141569, HD 169142, BD +40 4124, and LkHa224, see 
also lJuhasz et al.ll2010l ). The lack of the silicate emission feature suggests an unusually low 
density of submicron-sized silicate grains in the surface layer of the disk. One possible reason 
for this is that the small silicate grains settled to the midplane, do not r eceive any stellar 



irradiation, thus became too cold to emit at 10 fim (IMeeus et al.l l2002l : iDuUemond et al. 



20071 ). Alternatively, it may also happen that the silicate grains are too large to emit at 
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10 /iin due to grain growth (jCoUins et al.ll2009l ). 



Figs. [3] and show that the observed flux changes, if there are any, are within 30% 
and are mostly wavelength-independent: similar for continuum and PAH band wavelengths. 
In seeking a physical explanation for the observed variability, first we investigate whether 
the flux changes could be caused by fluctuating obscuration along the line-of-sight. Such 
a model is motivated by the fact that one of our targets, RRTau, is known to be UXor- 
type star. UXors, named after the prototype UXOri, show several magnitude deep optical 
minima due to variable circurn stellar obscuration along the line of sight in a nearly edge-on 
system (lEaton &: Herbstlll994j ). Theoretically, one could imagine that the observed mid-IR 
variability is caused by a "super-UXor" event when a circumstellar dust cloud obscures not 
only the star but the whole mid-IR emitting region. However, the observed magnitude of the 
mid-IR changes would require very large variations in the visual extinction (the extinction 
in the V band is ~30 times that at 5 /im), which is not observed at optical wavelengths for 
our targets. Moreover, the continuum and in particular the PAH emission originate from 
a very extended circumstellar region (at least several AUs in size), which is unlikely to be 
completely obscured. Therefore, we conclude that the UXor-phenomenon cannot be the 
physical origin of the mid-IR variations in PAH dominated objects. Any variability seen in 
the mid-IR regime should reflect variability in the emission from PAH or dust particles due 
to the changing illumination from the central source. 

Variable dust emission may occur when the stellar irradiation of the disk is modulated 
by time dependent shadowing effects. Changing the vertical structure of the inner disk can 
cast a shadow on a significant portion of the disk, decreasing its total stellar irradiation. 
However, this would only affect the continuum emitted by large grains and not the PAH 
emission. Since the PAH particles are strongly coupled to the gas and do not settle to the 
midpl ane, shadowing them is difficult, because the shadow should extend to a large sohd 
angle (IDullemond et al.l 120071 ). 



The only remaining way to change the irradiation of the disk is by varying the stellar 
luminosity. Herbig star s often show modest optical variability (see long- term optical moni- 
toring programs e.g. in lManfroid et al.lllQQll : iHerbst &: Shevchenkd 119991 ). whose amplitude 
is similar to the mid-IR changes we observed for the PAH dominated objects (Fig. H]). In 
some cases the optical variability is due to line of sight effects (extinction, long-lived stellar 
spots), which cannot change the luminosity of the star, consequently the integrated mid-IR 
emission would stay constant. In a few cases, however, it is possible that the central lu- 
minosity actually changes, for instance due to variable accretion. Herbig stars are modest 
accretors: accretion rates from Br7 observations are below a few times 10~^ Mq/jy. How- 
ever, there are indications that in some sources (e.g. RCrA), the Br7 equivalent width is 
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variable, pointing to variable accretion rate ( IGarcia Lopez et al.ll2006l ). The fact that not 



all intermediate mass stars exhibit luminosity changes is consistent with our finding that 
mid-IR variability is not widespread among PAH dominated stars (three variables out of 
eight). 

Our results suggest a circumstellar geometry where both the regions of PAH and contin- 
uum emission are directly illuminated by the central source. The most straightforward model 
to reproduce these results is a simple fiared disk geometry. In this model, the PAH particles 
and the large dust grains in the disk surface layer react to the changing illumination almost 
i nstantaneously. Even the optically thick disk midplane would react on a few years timescale 



(IChiang &: Goldreichlll997l ) , which is less than the temporal baseline of our variability survey. 
No opaque puffed-up inner rim is included in this model because no shadowing effects were 
revealed by our data. Note that our simple geometrical picture, dedu ced from the results 



of the variability measurements, is similar to the scheme proposed by iMeeus et al.l (120011 ) 
for their Group I sources of Herbig stars. The lack of small silicate grains high above the 
midplane of the disk, indicated by the weak or missing silicate emission feature, may help to 
keep clear the way for the stellar UV photons to reach the outer flared disk regions and excite 
the PAH molecules. In summary, varying illumination of a flared disk, possibly caused by 
variable accretion, is a promising idea to explain the mid-IR variability of PAH dominated 
objects. 



5.2. Objects with silicate emission 



The group of young stars exhibiting the 10 /xm silicate feature in emission contains both 
low- and intermediate mass objects: 20 Herbig Ae/Be and 16 TTauri stars belong here. 
Seven Herbig stars show PAH emission as well, but the PAH features do not dominate the 
spectrum (TTauri stars show no PAH emission in our sample). It is unclear why PAH 
emission is weak or absent in our intermediate mass objects with silicate emission. The 
distribution of spectral types among PAH dominated and silicate emission Herbig stars are 
similar, thus the reason might be related to e.g. intrinsically low PAH abundance in the 
disk atmosphere; the destruction of PAHs in the disk atmosphere by strong UV flux; o r 
shadowing or softening the stellar irradiation by the puffed-up inner rim (IKeller et al.ll2008l ). 
Taking into account that the PAH molecules are typically situated in the outer parts of a 
flared disk (Section l5.ip . the last p ossibility invokes a modestly flared or flat disk geometry, 
similar to the scheme proposed by iMeeus et al.l (120011 ) for their Group II Herbig stars. The 
fact that in Fig. [5] there is no obvious difference in the distribution of data points of silicate 
emission Herbig stars with and without PAH emission suggests that variability studies of 
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the silicate feature and continuum are not affected by the presence of PAH emission. This 
conclusion is supported by the flux ratios in Fig. |3] which suggest that, similarly to PAH 
dominated objects, for silicate emission Herbig stars the PAH features and the adjacent 
continuum show identical variability, thus magnitude changes in the 6-8 fim continuum and 
the 8-11.5 /im silicate feature domains are unaffected by the PAH features. Thus in the 
following we disregard the presence of PAH features in our silicate emission sources and 
concentrate on the 10 ^m silicate feature only. 



According to the standard picture of a two- layer passive disk (IChiang fc Goldreichlll997l ). 
the mid-IR continuum emission comes from warm dust grains both in the optically thick disk 
interior and in the superheated surface layer. The energy source is the central star, which 
illuminates the disk surface. In case of an actively accreting disk, there is an extra heating by 
the viscous energy released in the disk midplane, in the accretion columns, and in accretion 
hot spots on the stellar surface as well. The 10 /im silicate feature is thought to arise 
from the superheated silicate grains in the disk atmosphere. The radial extent of the disk 
atmosphere contributing to the mid-IR flux both in the continuum and in the 10 /im feature 
is within 10 AU for a low-mass star. The contribution of the disk interior is confined to 
within 1 AU. For intermediate mass stars, these dimensions are scaled up, an inner hole may 
be present, and the inn er rim of the disk, directly illuminated by the star, may be puffed-up 
(jPuUemond et al.ll200ll ). The rim may cast a s hadow on parts of the disk , strongly affecting 
the em itting area of the 10 /xm silicate feature. IVinkovic et al.l ( 120061 ) and lVinkovic fc Jurkic 
(120071 ) argue that many low-luminosity YSOs (both T Tauri and Herbig stars) have compact 
(^10 AU), optically thin, dusty halos around their inner disk regions. The dust temperature 
in the halo is set by the stellar heating. According to their modeling, the emission of the 
halo is dominant in the 2-5 /im regime, but it can also have important contribution at longer 
wavelengths. In certain cases (e.g. HP 163296) , there are claims tha t the halo alone suffices 
to explain the observed 10 /im silicate feature (IVinkovic et al.l 120061 ). 

The distributions of low- and intermediate mass silicate emission stars in Fig. [5] is 
somewhat different from each other. Most of the Herbig stars show very little variability 
in the 6— 8 /im continuum, while the 10 /im feature can vary as much as 0.6 mag (Fig. [5] 
upper right panel). On the other hand, T Tauri stars may vary also in the continuum 
(Fig. lower left panel). The number of T Tauri stars exhibiting wavelength independent 
flux change (diagonal strip) and showing larger changes in the silicate feature than in the 
continuum is comparable. The different variability patterns of low- and intermediate mass 
stars may indicate a difference also in their circumstellar structure. Another possible factor 
which might affect the variability statistics is the age of the objects: most low-mass stars in 
our sample are located in young star-forming regions, while the intermediate mass sample 
includes also isolated, somewhat older, and thus less variable, stars. 
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First we discuss possible physical reasons for the wavelength independent flux changes. 
We will basically consider the same scenarios as for the variability of PAH dominated ob- 
jects. Extinction of the mid-IR emitting region by a passing dust clump in the outer disk 
or envelope can be excluded because the accompanying optical flux changes would be un- 
reasonably high (up to 20 mag). Wavelength independent flux changes can be the result of 
varying illumination by the central star, possibly caused by time variable accretion onto the 
stellar surface. This scenario predicts similar amplitude changes also at shorter wavelengths. 
In a few cases, it was possible to test this by comparing two ISOPHOT-S spectra over the 
whole 2.5-11.6 fim wavelength range. Indeed, we found approximately wavelength indepen- 
dent behavior, for DRTau, Glass I, HD 104237, and EXLup. Although simultane ous optical 
photo metry does not exist, based on V-band light curves in the ASAS database (jPojmanski 
19971 ). the typical amplitude of the optical variability of these stars, which is often considered 
as a proxy for accretion variations, is similar to the observed < 0.6 mag mid-IR variability. 

For a surprisingly large number of our sources (including both TTauri and Herbig stars). 
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20081 : iMuzeroUe et al.ll2009l : iBary et al.ll2009l ). As we concluded in Section U]T1 the variation 
in the silicate emission is in all cases a change in the strength of the feature, and not a change 
in the spectral shape. 

For Herbig stars, an obvious explanation for the observed flux change s is that the inner 
rim c asts a shadow on parts of the disk. Using a very simple model, iDuUemond et al. 



( I2OOII ) investigated the effect of increasing rim height and found that it can cause increasing 
continuum and decre asing 10 iim feature , i.e. an anti-correlation between shorter and longer 
mid-IR wavelengths. lJuhasz et al.l ( 120071 ) performed a similar study using a radiative transfer 
code, and their model predicted a similar anti-correlation between 3—5 //m and >8 /im fluxes 
(although their target, SVCep, displayed a constant 10 //m flux, making it necessary to 
introduce an inner envelope in the system). ISitko et al.l ( 120081 ) proposed a slightly different 
model to explain mid-IR variability of the Herbig stars HD 31648 and HD 163296, suggesting 
structural changes in the gas disk close to or inside the dust sublimation zone. This gas 
disk may influence how much illumination reaches the inner rim from the central star. This 
model also predicts anti-correlation. Physical reasons for changes in the disk include thermal 
or magnetorotational instabilities and planetary perturbations. 



In certain cases, inner dust rims may also exist for TTauri stars (IMuzerolle et al.ll2003l ). 
and may play a role in the mid-IR v a riabil ity. Such a picture was discussed for LRLL31, 
a T Tauri star where iMuzeroUe et al.l (l2009l ) reported anti-correlation between shorter and 
longer mid-IR wavelengths, with a constant pivot point at 8.5 /xm. They proposed that, sim- 
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ilarly to Herbig stars, the height of the inner rim increases due to rising accretion luminosity, 
casting a shadow on the outer disk. Another scenario for LRLL31 is that the shadowing 
is due to dynamical perturbations by a stellar or planetary companion. These perturba- 
tions may induce war ps or spiral density waves in the inner disk, causing variable shadowing 
( iFlaherty et al.ll201ll ). CRCha (anti-correlation between A < 5 /xm and A > 8 //m), as well as 
WWCha and XXCha (anti-correlation between A = 6— 8/im and X>8fim) are promising 
candidates where one of these pher iomena might take pl ace. Note that spectral variability 
of CR Cha was also investigated by lEspaillat et al.l (120111 ) who found considerable change in 
the slope of the Spitzer/IRS spectra around 6 /im. 

The high variability amplitudes of the silicate feature of T Tauri stars might also be 
expl ained in the framework o f a recent model of turbulent disk accretion. [Turner et al.l (|2010[ ) 
and iHirose fc Turner! (120111 ) performed magnetohydrodynamical calculations of a stratified 
shearing box of the disk, and found that the magnetic activity intermittently lifts clouds 
of small grains into the disk atmosphere. The photosphere height changes by up to one- 
third over timescales of a few orbits, result ing also in changes of the mid-infrared surface 
brightness (see Fig. 18 in [Turner et al.ll2010l ). They also suggest that the changing shadows 
cast by the dust clouds on the outer disk are a cause of the daily to monthly mid-infrared 
variability found in many young stars. Since the dust clouds are optically thin, they may 
contribute primarily to the strength of the 10 /xm silicate feature, while they affect less the 
continuum emissi on. In the model, si gnificant variations in the 10 /im emission occur on the 
orbital timescale (ITurner et al.ll2010l ). In a typical T Tauri star, the disk area responsible 
for the bulk of the silicate emission is within the inner 1-2 AU, where the orbital period is 
less than a few years. Since our results on the variability also suggest timescales shorter 
than a few years (Section 14.21) . the Turner et al. model may be consistent with our findings 
concerning the variability timescale of the silicate feature. In order to check whether this 
mechanism could also account for the amplitudes of the observed v a riatio ns, we computed 
the predicted fluctuations in the following way. iHirose fc Turner! (120111 ) found that the 
starlight-absorbing surface moves up and down by a factor of two. The maximum possible 
10 /im variability amplitude will occur if unit starlight optical depth is reached where the 
temperature is near the 300 K needed to place the peak reprocessed emission at 10 /im. 
Under these favorable conditions, the 10 /im flux could change by a factor two, or 0.75 mag, 
similar to the mid-IR flux changes we observe in our silicate emission sample. We note 
that the disk photosphere height changes almost simultaneously at all azimuthal locations 
in the disk annulus in question, because the differential rotation rapidly shears out any 
non-axisymmetric structure, implying signiflcant changes in the emission of the system as 
a whole. Since the Turner et al. model predicts almost independent variability in the mid- 
infrared continuum and in the silicate feature, with higher amplitudes of the latter one, we 
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suggest it as a promising explanation for our observations. 



There might be other effects which change the sihcate feature and the adjacent contin- 
uum. One possibihty is that the system is binary or muhiple with muhiple disks. The sihcate 
features in each component may be different, and one of them may be variable, resulting 
in changes in the combined spectrum. Another option is what is happening for DG Tau, a 
single TTauri-type star: there is both em ission and absorption along the line-of-sight, and 
one of them is changing ( iBary et al.ll2009l ). There is also a possibility that the 10 fim silicate 
feature and the continuum changes independently from each other. This might occur if the 
continuum emission is dominated by viscous heating and the emission feature is dominated 
by illumination of the disk surface. The accretion in the disk midplane and the accretion onto 
the stellar surface changes independently, or at least with significant time difference. Finally, 
the example of CoKuTau/4 demonstrates that even for non- accreting sources, wavelength- 
dependent mid-IR variability can occur. This system consists of a close eccentric binary 
surrounded by a circumbinary disk. As the distance of the orbiting stars from the inner edge 
of the dust disk continuously changes in tir ne, the irradiation and consequently the mid-IR 
emission of the inner disk wall also varies (INagel et al.ll2010l ). Unfortunately, the present 
dataset does not allow us to decide between these scenarios. Future multi-epoch simulta- 
neous optical/near-IR/mid-IR observations, and detailed modeling of each source would be 
the key to conclusive results. 



5.3. Objects with silicate/ice absorption 

Stars with silicate and occasionally with ice absorption bands are all low-mass stars in 
our sample. The increasing SED and the absorption features indicate that these objects 
are deeply embedded in a cold, dense, extended envelope. The mid-IR emission is probably 
coming from a circumstellar disk as well as from the inner warm part of the envelope. The fact 
that they are mostly distributed along the diagonal stripe in Fig. [5] suggests scenarios that 
cause wavelength-independent flux changes. The short wavelength (A < 5 /xm) flux changes, 
where available, support this idea. 

In accordance with this result, the concept of shadowing as discussed for disks is not 
applicable for envelopes. This is not surprising, since the envelope covers a large solid angle, 
thus the shadowing material within the inner edge of the envelope should cover a large solid 
angle. In principle, shadowing effects might occur in the disk, but these changes would 



be te mpered by the backwarming of the outer disk by the surrounding envelope ( iNatta 



19931 ). Similarly to PAH dominated and silicate emission objects in the diagonal stripe, for 



silicate absorption objects, the most obvious explanation for the varying central luminosity is 
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variable accretion rate. A good example for this is the well-studied triple system T Tau, where 
the Sa com ponent is responsible for the 10 /xm silicate absorption and for the mid-infrared 
variability. Ivan Boekel et al.l ( 120101 ) show that the short-term (daily-weekly) mid-infrared 
variability of T Tau Sa can be attributed to variable accretion, while the long-term (yearly) 
flux changes can be either due to variable accretion, variable foreground extinction, or the 
combination of these two mechanisms. 



6. Summary and conclusions 

In this paper we present low-resolution mid-IR spectra of 68 low- and intermediate 
mass young stars obtained with ISOPHOT-S between 1996 and 1998 and with Spitzer/IRS 
between 2004 and 2007. Utilizing self-developed software packages, we interactively re- 
processed the spectra for improved removal of instrumental artifacts. By comparing multi- 
epoch spectra of each object, we analyzed the statistics and the trends of mid-IR spectral 
variability and its implications on the geometry of the circumstellar material. Our study has 
three novel aspects: (1) the spectroscopic observations enable the study of spectral variations 
in a relatively extended mid-IR wavelength range; (2) instead of focusing on a certain star 
forming region, individual well-known YSOs were observed; (3) and the timescale of our 
ISO/Spitzer comparison is longer than previous mid-IR variability studies, covering a full 
decade. Our main achievements are the following: 

• The mid-infrared spectral atlas constructed for this project presents an improved, final 
data reduction of ISOPHOT-S spectra. The atlas contains all low- and intermediate 
mass stars ever observed with the ISOPHOT-S instrument, thus, it constitutes the 
legacy of the ISO. 

• Our spectra show that mid-IR variability among low- and intermediate mass YSOs 
is ubiquitous. We calculated synthetic photometry in the 6-8 //m and the 8-11.5 ^m 
wavelength range, and found that 79% of the sources vary more than 0.1 mag, while 
43% is variable above the 0.3 mag level. A comparison of the variability characteristics 
on annual and decadal timescales revealed no significant differences, constraining the 
physical timescale of the variability to less than a few years. 

• For intermediate mass stars with spectra dominated by PAH emission, we found rela- 
tively low incidence of variability. Flux changes in this group are mostly wavelength- 
independent, and can be interpreted in terms of non-steady irradiation of the disk due 
to fluctuating accretion. We propose a simple flared disk geometry to model these 
sources. 
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• Intermediate mass stars exhibiting silicate emission at 10 /im often show higher vari- 
abihty amphtude in the sihcate feature than in the adjacent continuum. Shadowing 
the disk by vertical variations of a puffed-up inner disk rim is invoked to explain the 
wavelength-dependence of the flux changes. The deduced geometry is a modestly flared 
or flat disk. 

• TTauri stars are the most frequent variables in our sample. For those which exhibit 
more pronounced variability in the 10 /xm silicate feature than in the continuum we 
propose to consider a new model: strong magnetohydrodynamical turbulence may 
intermittently lift clouds of small gr ains into the disk atmosphere resulting in extra 



silicate emission (ITurner et al.ll2010l ). Nevertheless, based on our spectra alone, we 



cannot exclude that other scenarios may also be applicable. 

Sources exhibiting silicate or ice absorption are objects embedded in a dense envelope. 
They typically show wavelength independent flux changes, probably due to varying 
accretion rate. 



Our results suggest that mid-IR variability is widespread among YSOs. Interpreting the am- 
plitude, wavelength-dependence, and timescales of these flux changes is a new and promising 
possibility to explore the structure of the circumstellar disks and their dynamical processes. 
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A. Interactive post-processing of ISO observations 

Though the ISO Legacy Archive contains rehable results for most ISOPHOT-S observa- 
tions, several instrumental artifacts were discovered after the closure of the Archive in 2001. 
We have developed an IDL-based processing package aiming at correcting these artifacts. 
Possible correction algorithms were tested and optimized on a set of 43 normal star mea- 
surements, and the same data set was used to estimate typical measurement uncertainties 
via comparison with photospheric models. The processing scheme consists of the following 
steps: 



Deglitching. As a preparation, we checked whether PIA had successfully removed all 
cosmic glitches from the data. Although the PIA built-in deglitching algorithms are quite 
efficient, we often experience that glitches are still present in observations of very faint 
sources, producing artificial spikes in the spectra. The reason is that at low signal level the 
recovery from a cosmic hit takes a relatively long time, and the disturbance following the 
hit may have a characteristic timescale comparable to the full measurement time (for an 
example see Fig.E]). In such a situation, the statistical algorithms of PIA - which assume 
that most of the data points are not affected by glitches - may not work properly. In order 
to identify the remaining glitches, the signal evolution of each pixel was visually inspected at 
the intermediate Signal Raw Data (SRD) processing level, and - if necessary - the affected 
data points were manually discarded using PIA. 



Detector temperature. It was suggested during the ISO mission that when the detector 
temperature was outside the nominal 2.8-3.1 K range, the transient behavior might have dif- 
fered from the normal one, leading to systematic photometric errors ( lAcosta-Pulido fc Abraham 



20031 ) . We checked all observations presented here, and found that one measurement (HD 142666, 



TDT: 10402847) had non-nominal detector temperature. In the case of this measurement 
we computed for every pixel the dispersion of the calibrated flux values derived from the 
individual observing ramps by the dynamic calibration procedure. We checked whether the 
dispersion was higher in this affected measurement than in other data sets of similarly bright 
stars measured in the nominal temperature range. We indeed found a 2-3cr increase in the 
noise, however, it is only 3-4% and 1-2% of the total measured flux in the ISOPHOT-SS 
and ISOPHOT-SL wavelength range, respectively. Thus, we conclude that non-nominal de- 
tector temperature is not a major source of uncertainty in the ISOPHOT-S observation of 
HD 142666, and we will neglect this effect in the processing. 



- 25 - 



Orbital phase. It was suspected that observations obtained either at the very beginning 
(orbital phase < 0.25) or at the very end of the orbit (orbital phase > 0.80) might have 
suffered from reduced accuracy due to a higher dark current and to a higher cosmic glitch rate, 
respectively. In our sample, 9 measurements were taken during early, and 18 measurement 
during late orbital phase. In order to check if there is an increase of uncertainty related to the 
dark current subtraction at early orbital phases, we tested the distribution of dedicated dark 
current calibration measurements around a fitted curve that represents the orbital variation 
of the da rk current per pixel, and found n o higher scatter at early orbital phases (see also 
Fig. 1 in lAcosta-Pulido fc Abrahamll2003l ). This means that the subtraction of the orbital 
dependent dark signal does not introduce a higher noise component at early phases than 
later in the orbit. Close to the end of the science window of the orbit, the rate of cosmic 
glitches could somewhat rise, possibly making the data noisier than at earlier orbital phases. 
However, this noise is probably random in nature, and there is no reason to assume correlated 
behavior of groups of pixels. Since our analysis is based on synthetic photometry derived 
from the spectra, our results are unaffected by the increased noise at late orbital phase. 



Signal memory from preceding observation. We found the appearance of artificial 
spectral features in some spectra due to memory effects from a preceding observation of a 
bright source (Fig. [7]). In order to identify affected spectra, we studied the short (32 s) dark 
measurements performed before each ISOPHOT-S observation. After reducing the complete 
sample of these measurements from the mission, we determined their average signal levels 
and typical measurement uncertainties per pixel (dashed and dotted lines in Fig. |8]). If the 
measured short dark signal of a particular observation exceeded the average level by more 
than la for a continuous section of the spectrum (i.e. for a group of neighboring pixels, for 
an example, see Fig. |8]), then the observation was flagged for memory effect. 

From tests on cases where external information on the spectral shape (theoretical mod- 
els, TIMMI20 observations) was available, we concluded that subtracting the excess of the 
short dark signal from the astronomical observation per pixel can efficiently correct for the 
memory effect. In order to determine the accuracy of this correction, we analyzed the only 
standard star observation that was strongly affected by signal memory effect. We found that 
for pixels where the correction was below 50% of the measured signal, the uncertainty of the 
corrected signal remained below 8-10%, i.e less than the photometric uncertainty assigned 
to the measurements in the final step of our processing scheme. In the 50-100% memory 
excess range, the precision of the corrected fiuxes was 20%, and even pixels suffering from 



t Thermal Infrared Multimode Instrument 



- 26 - 



very high memory contamination, in the range of 100-300%, could be corrected to have a 
standard deviation of 33%. In our sample, we performed memory correction on 19 observa- 
tions of 17 different sources (marked with M in Column (7) of Tabled]). In almost all cases 
the correction was in the 5-20% range, and the associated uncertainly was probably lower 
than the general photometric accuracy. The two more strongly affected cases (marked with 
M! in Tabled]) are CedllllRSS (TDT 62501412), and CK2 (TDT 10803228). Out of these, 
the SL part of the spectra of Ced 111 IRS 5 could be used for variability studies, while CK 2 
was discarded from our analysis. 



Background subtraction. For 61 observations out of the 94 ISOPHOT-S spectra pre- 
sented here, no corresponding sky background measurement has been performed. In the 
wavelength range of ISOPHOT-S the dominant background component is the zodiacal light 
(apart from localized regions like HII regions, reflection nebulae, or the galacti c plane). The 



zodia cal flux was clearly measurable for ISOPHOT-S in the 5.8-11.7 fim regime (ILeinert et al 



20021 ) . and its contribution is not negligible for fainter stars (Fig.lH]). We created a model 
which was able to predict the spectrum of the zodiacal background towards a given direc- 
tion and on a given date. The algorithm used 4.9 and 12 //m photometric points from 
COBE/DIRBE photometry as input, fitted a Planck curve to them, and finally scaled 
with a wavelength- and exposure time-dependent factor optimized for t he 29 high qual- 



ity ISOPHOT-S spectra of the zodiacal light presented in iLeinert et al.l (120021 ). Those 61 
measurements where background was predicted in this way are marked with B in Column 
(7) of Table dl 

We also performed an error analysis of the correction by computing a predicted back- 
ground spectrum for those 18 staring observations where a background measurement was 
available, and compare the two spectra pixel-by-pixel. The results demonstrated that the 
average error bar for the ISOPHOT-SS part is between 0.04 and 0.07 Jy, while for the 
ISOPHOT-SL domain is between 0.03 and 0.13 Jy, monotonically increasing towards longer 
wavelengths. For most of our targets, this additional uncertainty is well below the over- 
all 10% precision of the ISOPHOT-S calibration. For the five faintest sources (among those 
where predicted background was subtracted), however, the uncertainty of the predicted back- 
ground exceeds 10% of the measured flux. These cases are BARN 5 IRS 1, LDN1551IRS5, 
Ced 111 IRS 5, [SVS76]Ser4, and CK2. For the first three objects, only the ISOPHOT-SS 
part is affected, while for the latter two, the whole ISOPHOT-S spectrum is uncertain. In 
Column 7 of Table dl we mark six measurements of these five objects with an exclamation 
mark (B!). We conclude that variability can still be discussed for the first three objects 
using the ISOPHOT-SL part (and Spitzer when available), while for the two latter cases, 
variability study is not possible. 
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Special processing of spectral maps. There are a few measurements which were per- 
formed in mapping mode (usually ID scans). The target is usually situated on the central 
map position, and the neighboring raster steps can be used to obtain a background spectrum 
interpolated for the source's position. The mapping mode has several advantages over the 
typical ON-OFF staring observations, because it provides a well-measured background, and 
the subtraction of the interpolated background automatically cancels any memory effect, 
thus the algorithm described above is not used here. It is a problem, however, that the 
dynamic response calibration method (lAcosta-Pulido fc Abraham! |2003| ) was developed for 
staring observations, and maps are standardly calibrated by means of a static spectral re- 
sponse function. This method, however, cannot cope with the different transient timescales 
of the different detector pixels, and may produce obvious spectral artifacts. In order to 
find an alternative solution, we tested a simple algorithm which considers the time sequence 
of the map as a long staring observation, applies Dynamic response calibration on it, and 
creates a map from the calibrated flux values. The results were in good agreement with 
the ISOPHOT photometry and lacked any obvious spectral artifacts. Raster observations of 
LDN 1489 IRS and 00 Ser were processed using this algorithm. 



Off-center position of the source in the beam. The footprints of the ISOPHOT-S 
pixels were peaked rather than fiat-top, therefore observing a compact source outside the 
optical axis of the pixel could change the measured signal. Since the calibration was set up 
for the pixel center, off-center observations have to be corrected in order to derive correct 
flux values. In addition, because the footprint profiles varied with wavelength (there were 
especially large jumps at 3.7 and 8.8 fim) off-center position led to the appearance of spectral 
artifacts (an example is shown in Fig. fTUl) . In order to correct the flux values and the spectral 
artifacts, we determined the offset of the source from the center by comparing the accurate 
position of the objec t with ISO's coordin ates. The 2-dimensional footprint map (we utilized 



the ones derived by iLeinert et al.l |2002| ) was then sampled at the offset position, and the 
measured flux was scaled with the ratio between footprint values at the center and at the 
offset location. This procedure was repeated for each detector pixel independently. 

In order to estimate the uncertainty related to the correction, we derived error bars for 
each point of the footprint map by comparing those two original footprint measurements 
of the calibration star HR 7924 which were averaged to create the final footprint map. We 
found that in the majority of the measurements the estimated error was less than 5% of the 
total flux in the whole ISOPHOT-S wavelength range. Offset correction was applied to 66 
ISOPHOT-S spectra of 50 different sources (marked with O in Column (7) of Tabled]). There 
are only two seriously affected (> 10%) data sets: the short wavelength part of one of the 
two spectra of WW Vul (TDT: 17600465) and the spectrum of MWC863 (TDT: 28900460). 
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Both measurements are flagged with an exclamation mark (O!) in Table [H and are discarded 
from the variability analysis. 

Empirical photometric correction. In order to check for any remaining systematic 
effects in the final photometry, we queried the ISO Archive, collected all normal star obser- 
vations, and reduced them in the way described above. Fig. [TT] shows a typical result, where 
the measured spectrum is compared with the expected one. The latter spectrum was taken 
from the ISO calibration data base which includes a compilation of predicted spectra pro- 
vided by M. Cohen or P. Hammersley (http:/ /iso.esac.esa.int/users/expLlib/ISO/wwwcal/). 
For those objects not included in the data base we took one of the available models of a star 
of identical spectral type, and scaled the values to the K-band magnitude of our object. In 
our sample of 43 observations of 28 normal stars, no excess due to hot circumstellar dust is 
expected (the sample also includes the calibrator stars). 

In Fig. [12] we show the measured-to-predicted ratios of 3 representative pixels for the 
whole normal star sample. For each pixel, the ratios above a brightness threshold (usually 
5 Jy) were averaged and their standard deviation was computed. For the faint stars (below 
1 Jy), the ratio values became too noisy, and we computed the average and the standard 
deviation of differences between the measured and predicted fluxes. The averages derived 
this way represent the typical systematic errors. We decided to correct for these systematic 
offsets by dividing by the average flux ratio (for bright sources) or by subtracting the average 
flux difference (for faint sources) from all observations. This empirical photometric correction 
was the last step of our post-processing sequence. 

Error budget. The formal uncertainties provided by the PIA reduction do not contain 
several important sources of uncertainties, e.g. the ones related to the corrections in the 
post-processing. The average uncertainties in Fig. [131 the other hand, do not represent 
the quality differences among observations (e.g. the effect of a particularly energetic cosmic 
glitch). As final photometric uncertainties, we decided to adopt the maximum of the two 
types of error bars. In practice, for most pixels/measurements the average uncertainty 
values are the higher. Their typical values are 5-10% for bright sources, 30-100 mJy for 
faint sources. 
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B. Notes on individual sources 



In the following, we collected literature information on some of the sources, in order 
to help the reader with the understanding and the interpretation of the spectra plotted in 
Fig. |2j We mostly focused on three points: information on (1) variability in the optical- 
IR wavelength regime, (2) multiplicity, and (3) extended emission at optical-IR wavelengths. 
Single sources with no extended emission, and whose observations are unaffected by technical 
artifacts are not discussed here. Note that here we also included sources where variability 
could not be investigated. 



IRAS 03260+3111 (also known as NGC 133 3 SVS 3) is a young Herbig binary illuminat- 
ing a small reflection nebula (IStrom et al.lll976l). T he binary is composed of a B5 and an F2 
star, with a separation of 3'.'62 ( IStraizys et al.ll2002l ). The s econdary is 32 time s fainter than 



the primary in the K band and 24 times in the L band. ( Haisch et al. 



appears extended with a fan-shaped nebulosity at 10 /im (IHaisch et al. 



2004 ). The system 
20061 ). There is no 



evidence of emission from silicate grains in the mid-IR spectrum of IRAS 03260+3111. On 
the other hand, PAHs emission is observed not only at the position of the Herbig star, but 
also at various locations in the NGC 1333 nebula. The different PAH features change with 
changing distance fr om the central source due to the changing relative populatio r is of ionized 
and neutral PAHs ( iBregman et al.lll993l : iRoche et al.lll994j : iJoblin et al.lll996l : ISloan et al. 



19991 ). The MIR spectra published in these papers all have different flux levels, and they 



are also different from our ISOPHOT-S spectrum. These differences are due to the different 
beam sizes of the instruments and the exact pointing. 



BARN 5 IRS 1 is a low-mass embedded young stellar object in the Barnard 5 dark cloud 
( iBeichman et al.lll984j]. It is surrounded by a v ery faint reflection nebula detectable at near- 
infrared wavelengths ( iMoore fc EmersonI 119921 ). We found no evidence in the literature for 
extended mid-infrared emission. JHKL ph otometric monitoring revealed that the source has 
faded significantly between 1983 and 1993 ( iMoore fc Emersorull992l . Il994j ). In the case of the 
ISOPHOT-S spectrum, predicted background was subtracted, making the short wavelength 
part more uncertain. However, the SED of the source rises steeply towards longer wave- 
lengths, and the long-wavelength part of the ISOPHOT-S spectrum can be safely compared 
with the Spitzer/IRS spectrum. 



LDN 1489 IRS is a single embedded YSO ( IConnelley et al.ll2008[ ). The source looks ex- 
tended in the 2MASS J, H, and K images, as well as in the 3.6 fim Spitzer/IRAC image. 
Although it is point-like in the IRAC images at >5.8 //m, we cannot exclude the possibility of 
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extended emission, which may be the reason for the shght difference between the ISOPHOT- 
S and the Spitzer/IRS flux levels. Thus, this source is not included in the discussion about 
variability. 



T Tau is a triple system consisting o f T Tau N, T Tau S a, and T Tau Sb. The largest sepa- 
ration between the components is 0'.'7 (IKohler et al.ll2008l ). therefore neither ISO nor Spitzer 
resolves the system: the observed flux represents the sum of all c omponents. T TauS a 
is a Herbig star, wh ile TTauSb and TTauN are low- mass objects (IDuchene et al.ll2006l ). 
Skemer et al.l (120081) presented photom etry at 8.7, 10.55, and 11.86 fim for each of the three 



components, while iRatzka et al.l (120091 ) show 8-13 /im spectra for T Tau S and T Tau N. Based 
on these results, it seems that only T Tau Sa displays silicate absorption, while T Tau N and 
TTauSb h a ve mo re-or-less featureless 10 /im spectra (TTauN may display slight emission). 
Beck et al.l (j200J) resolve T Tau N and T Tau S in the 2-4 /im range, and claim that T Tau S 
shows deep, broad H2O ice absorption at 3.1 /im, while TTauN is featureles s in th is wave- 
length range. Photome tric monitoring between 2 and 12 fim by iBeck et al.l (120041 ) and by 
van Boekel et al.l (120101 ) indicate that TTauN is constant within the uncertainties, while 
TTauS shows significant variability. 



DGTau 



IS a sme 



ingle classical T Tauri star ( Leinert et al.lEoQl : Connelley et al. I l2008h . IWooden et al 



(I2OOOI ) report on multi-epoch mid-IR spectra where they found that the 10 /im silicate feature 
showed significant variability on time scales of months to years, while the underlying contin- 
uum remained relatively constant. They claim that prior to 1996, DGTau had a featureless 
10 fim spectr um, which then tu rned into silicate emission, then went into absorption at the 
end of 1998. iBary et al.l (120091 ) also observed that the silicate feature displayed variabihty 
on monthly and yearly timescales. All our ISOPHOT-S and Spitzer spectra are featureless, 
but they indicate that the continuum is variable. 



Haro6-10 (also kno wn as GVTau) is a T Tauri binary with a separation of 1'.'2, P.A. 355° 
(ILeinert fc Haaslll989l ). Both the ISO and the Spitzer spectra contain the full fiux from both 
components (the Spitzer slit is oriented along the binary). At wavelengths shorter than 4 /im, 
the southern component dominates , while most of the em i ssion above 4 / im corn es from the 
extremely red northern component (ILeinert fc Haaslll989l ). iLeinert et al.l (l2001al ) reported J, 
H, K, L', and M-band variability, as well as the variation of the 3.1 /im ice-band absorption. 
Their resolved observations prove that both components vary significantly on timescales of 
a few months. Our ISO and Spitzer spectra confirm that the system is variable. 
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LDN 1551 IRS 5 is a deeply embedded low-mass triple system consisting of a binary with 
projecte d separation of (y/3, and a third component at a distance of (/.'OQ of the northern com- 
ponent ( iLim &: Takakuwall2006[ ). each star harboring a small circumstellar dis k. The whole 
syste m is surrounded b y a circumbinary disk and a flattened envelope as well ( lOsorio et al. 
20031 ). iLiu et al.l (119961 ) observed the source at 12 /im, and by comparing the observed flux 
with previous 12 and N-band photometry, reported significant mid-IR variability. Al- 
though LDN 1551 IRS 5 displays some extended emission at all Spitzer/IRAC images, the 
larger aperture ISO spectrum has a lower flux than the Spitzer spectrum, indicating real 
variability. We note that in the case of the ISOPHOT-S spectrum, predicted background 
was subtracted. The flux of LDN 1551 IRS 5 steeply rises towards longer wavelengths, but be- 
low 5 ^m, the source is so faint that the uncertainty of the background subtraction dominates 
the error budget, and the overall error is above 10%. 



HL Tau is a single classical T Tauri star. At optical wavelengths, the object is a compact 
reflectio n nebula with complic ated morphology and d iameter of ~3^^; no po int source is 
visible (IStapelfeldt et al.l 119951 ). In the JHK images of iMurakawa et al.l (120081 ). the central 
source can be seen, but it is still surrounded by extended nebulosity (about 3" in extension). 
According to our knowledge, there is no information on the literature on whether the source is 
extended at mid-IR wavelengths. The fact that the larger aperture ISO spectrum has a lower 
flux level than the Spitzer sp ectrum indicate rea l varia bility. The 8-13 yum UKIRT/CGS3 
spectrum obtained in 1993 by lBowey &: AdamsonI (120011 ) is consistent with the ISOPHOT-S 
brightness level. 



UZ Tau is a hierarchical quadruple T Tauri system. UZ Tau E is a spectroscopic binary 
with an orbital p eriod of 19.1 days. At a distance of 3'.'8 can be found UZ Tau W, which itself 
is a 0'.'34 binary ( iPrato et al.ll2002l . and references therein). Due to periodic accretion from 
the circumbinary disk, UZ Tau E shows optical variability o f about 0.4-1.0 mag in the V, R, 



and I bands, while UZTauW is constant within 0.02 mag ( 1 Jensen et al.l 120071 : iKospal et al. 



201ll ). UZTauE is also variable at millimeter wavelengths, possibly related to non-thermal 
emission from magnetosph eric reconnection events above the binary components' photo- 
spheres (IKospal et al.l 1201 ll ). The Spitzer/IRS slit was centered on UZTauE, although the 
W component also has some contribution to the observed flux. The ISOPHOT-S measure- 
ments were centered between the E and W components, and essentially measured the sum of 
the two components' fluxes. For this reason, the Spitzer/IRS and ISOPHOT-S spectra can- 
not be compared, but variability can be studied by comparing the two ISOPHOT-S spectra, 
which indicate some flux changes in the silicate feature around 10 fim. 
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VY Tau is a subarcsecond binary with a separation o f (y/66, P.A. of 317°, K - band flux ratio 



of 3.8 , and L-band flux ratio of 6.1 (ILeinert et al.lll993l : iMcCabe et al.ll2006[ ). iRichichi et al. 



( Il994j ) note that the K-band brightness ratio of the two stars is time- variable. VYTau is 
an EXor-type object, exhibiting sporad ic l-4niag optical outbursts, but has b een in quies- 
cence at V~13.8mag since about 1972 ( lHerbiglll977 . 



1990 



Herbst et al. 1994). The source 



shows photospheric emission in the JHK bands ( IShiba et al. 



1993 



Sipos et al.l 1201 ll ). but 



Furlan et al. 


2006; 


Sargent et al. 


2009; 


Watson et al. 


2009; 


Sipos et al. 


2011) 



indication in the literature for extended IR emission. Both the ISO and the Spitzer spectra 
contain the full flux from the binary components, thus they can be compared. 



DR Tau is a highly accreting single classical T Tauri star with a long history of optical and 



near-IR photometric and spectrosc opic variab i 



with fluctuating accretion rate fe.g. iGotz 



1980; 



ity, which point to magnetospheric accretion 



Isobe et al 



1988; Hessman fc Guenther 1997; 



Smith et al.lll999l : lAlencar et al.ll200l[ ). iKenyon et al.l (Il994[ ) present light curves from the B 



to the L band, and explain the observed flux variations with the presence of a hot spot (the 
base of the accretion flow) on the stellar surface. PR Tau is sometirnes cla ssifled as an EXor, 
a young star exhibiting short accretion outbursts (ILorenzetti et al.l 120091 ). Our observations 
indicate that DRTau is also variable in the whole 2.5-15 //m wavelength range. 



HD 34700 is a quadruple T Tauri system. The primary itself is a double-line spectro- 
scopic binary (HP 34700 Aa and HP 34700 Ab ), consisting of two equal- mass GO- type stars 
f lArellano Ferro fc Giridhai]l2003[ lTorresll2004l ) . Two faint stars at a distance of 5'.'2 (HP 34700 B, 
spectral type: Ml-2, classical T Tauri) and 9'.'2 ( HP 34700 C, spectr al type: M3-4, weak-hne 
T Tauri) are also associated with HP 34700 A ( ISterzik et al.ll2005l ). The IRS slit only in- 
cludes the spectroscopic binary HP 34700 A, while the ISOPHOT-S aperture includes all 
four components. Although HP 34700 B seems to be actively accreting, due to their faint- 
ness and distance from component A, components B and C have negligible contribution to 
the ISOPHOT-S flux. Thus, both instruments essentially measure HP 34700 A. 



RRTau is a Herbig Ae star. iGrinin et al.l ( 12002| ) observed a deep optical and near-IR 
minimum in the light curve of RR Tau in 2000-2001. They found that the UBVRIJH fluxes 
showed synchronous behavior, and the color changes indicated that the dimming was caused 
by a circumstellar dust cloud crossing the line of sight. However, the fading in the optical and 
JH bands were accompanied with a flux increase in the K and L band, suggesting that flux 
changes at longer wavelength have a different physical mechanism than variable extinction. 
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Reipurth SON IRS 1 . Reipur th 50 was first reported as a nebulous object by iReipurth 
( 119851 ) ■ IReipurth &: Ballyl (Il986[ ) reported the sudden appearance of a highly variable, coni- 



cal ne bula Reipurth 50 N, possibly due to an FU Ori-like eruption. IScarrott &: Wolstencroft 



( 119881 ) identified Reipurth 50 N IRS 1 (a.k.a. HBC494) as the illuminating source of the 
Reipurth 50 N refie c tion nebula. Near-IR imaging, imaging polarimetry, and CVF spec- 
troscopy by ICasalil ( 1l99ll ) revealed that there is only one single source in the Reipurth 50 N 
cloud, IRS 1, which is a point source at 3.6 /xm. 



SX C ha is a T Tauri-type binary with a separation of 2'.'2, and P.A. of 288° ( iReipurth fc Zinnecker 
1993h. The east and w est components have spectral types of MO. 5 and M3.5, respectively 
( jMuzeroUe et al. 2005a). The 2MASS images show that while the west component is brighter 
in the J band, the east component is brighter in the K band. The Spitzer slit is centered on 
the east component, and the west component is 1'.'2-1'.'6 away from the slit centerline. This 
means that 55-70% of the fiux of the west component is also included in the Spitzer measure- 
ment. The ISOPHOT-S data, however, contains the total fiux of both components. Thus, 
due to beam confusion, the spectra of the two instruments cannot be directly compared. 



HP 9 5881 is a Herbig Ae star. Based on spectro-astrometric observations, iBaines et al. 



( 120061 ) consider it a possible subarcsecond binary. We found no indication in the literature 
for extended emission at any wavelength, thus the ISOPHOT-S and Spitzer spectra can be 
safely compared. The two spectra agree within the measurement uncertainties until about 
10 /xm. Between 10 and 11.5 /im, the ISOPHO T-S spectrum is h i gher t han the Spitzer at a 
la level. The TIMMI2 spectrum published by Ivan Boekel et al.l (|2005[ ) is even higher than 
the ISOPHOT-S, indicating possible variability. 



CT C ha is a single T Tauri-type star (iGuenther et al.l 120071 ). Recently, ISchmidt et al. 



( 120081 ) found a candidate sub-stellar companion (brown dwarf or planet) by direct imag- 
ing, at a projected separation of 2'.'67. Since the companion is very faint (K-band flux ratio 
is 313), it does not affect our spectra, and the ISOPHOT-S and Spitzer/IRS data can be 
safely compared. 



VW Cha is a triple system. The separation of the A and B components is 0'.'66, P.A. 
177°, and their K-band flux ratio is 25.6. The se paration of the B and C components is 
O'.'IO, P.A. 233°, and their K-band flux ratio is 1.34 ( iBrandeker et al.ll200ll). There is also a n 
even fainter, fourth component at a distance of 16.78" from VWChaA (ICorreia et al.ll2006l ). 
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Both our ISO and Spitzer spectra contain the flux from the three brightest components (A, 
B, and C). 



CU Cha (also known as HP 97048) is a single Herbig star with spectral type of B9.5/A0 



Ghez et al.lll997l : lBailevlll998l : lDoucet et al.ll2nn7l ). Based on ISOCAM-CVF data. ISiebenmorgen et al 



: bund that the emission of CU Cha in the PAH bands is extended on scales of about 



5-10". According to iDoucet et al.l ( 120071 ). this emission comes from transient l y hea ted PAH 
molecules in an extended envelope around the star-disk system. iLagage et al.l ( l2006l ) claimed 
that CUCha is extended in the 8.6 /jm PAH band, and the brightness profile is well repro- 
duced by emission from the surface layer of a vertically opti c ally t hick, inclined, flared 



disk, extending at least up to 370 AU (or 2'.'1). iDoucet et al.l (120071) found that the ob- 
ject is also extended in the 11.3 //m PAH band, while Ivan Boekel et al.l (120041 ) could resolve 
the disk also at continuum wavelengths. While the differences between the ISOPHOT-S 
and Spitzer/IRS spectra might be due to the different beam sizes, the fact that the three 
ISOPHOT-S spectra differ by 20-40% points to real variability. Long-term op tical monitoring 



indicates that the amplitude of op tical variability of CU Cha is < 0.1-0.3 mag (IManfroid et al. 
199ll : iHerbst fc Shevchenk3ll999l ). 



Glass I is a T Tauri-type binary with a separation of 2.67", P.A. of 285°, and K-band 
flux ratio of 5.9 (IHaisch et al.l 120041 ). It is also known as CedlllIRS4 or ChaIT33. The 
primary is brighter in the optical, but it is a "na ked" TTauri; most of the IR emission come s 
from the deeply ern bedded Class I secondary (iFeigelson &: Krisd Il989l : iPrusti et al.lll99ll ). 
Giirtler et al.l (Il999[ ) already reported strong mid-IR variability based on two ISO spectra. 
The Spitzer slit in 2006 was positioned halfway between the two stars, but the slit and the 
P.A. of the binary is misaligned, leading to some flux loss. We do not attempt to correct 
for this effect. The slit in 2008 was positioned exactly on the secondary component, which 
dominates the IR emission. Our variability study for this source is based on the ISOPHOT-S 
spectra, which contain the total flux of both components. 



CedlllIRS5 is a single TTauri-type star (IHaisch et al.ll2004j ). It is extended on an 
arcminute scale, even at the 8 /zm IRAC image. It is also called Chamaeleon IR nebula. 
This explains the difference between the ISO and Spitzer spectrum. However, variability 
can be studied by comparing the long-wavelength parts of the two ISOPHOT-S spectra, 
which have nearly identical pointings and are not strongly affected by memory effects or 
uncertainties related to the subtraction of the predicted background. 
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HD 97300 is a B9-type binary Herbig star ( ISiebenmorgen et al.lll998l ). In the HR diagram. 



it is l ocated close to the ZAMS, and it has no significant IR excess below 24 fim ( iKospal et al. 



2012h. The separation of the binary co n aponents is 0^/8, and their K-band fiux ratio is 17 



( IGhez et al.lll997l : iLafreniere et al.ll2008l ). ISiebenmorgen et al. Jl998h noticed that at certain 
IR wavelengths, two emission peaks can be seen, one coinciding with HD 97300, one being at 
a distance of about 3". They also detected an extended, ring-like structure, whose spectrum 
is dominated by PAH emission. The size of the ring is about 50"x33". It is not centered 
on the star, and is probably made of interstellar matter, whose density is enh anced by the 



interaction of the stellar wind from HD 97300 with the environment (see also iKospal et al. 



20121 ). The PAH molecules in the ring are transiently heated by HD 97300. This extended 



emission explains the factor of 3-4 difference we de tect between the ISOPHO T-S and IRS 
spectra (the ISOCAM-CVF spectrum published by ISiebenmorgen et al.l l2000l . representing 
the total flux of the system, is another factor of 2 higher than the ISOPHOT-S spectrum). 
The good agreement between the two ISOPHOT-S spectra indicates that the emission of the 
central source, and consequently that of the PAH molecules heated by it, is constant in time. 
Indeed, long-term optical and near-IR photo metric monitoring indicates that the brig htness 
of HD 97300 is constant within 0.1-0.2 mag jManfroid etal]ll99ll : bavies et al.lll99o[ l. 



Cedll2IRS4 (also known as FMCha) is a single TTauri-ty pe star with no indicatio n 
of extended near - IR or mid-IR emission in its immediate vicinity ( iHaisch et al.ll2004l |2006| ). 
Alexander et al.l ( 120031 ) presented an ISOCAM-CVF spectrum of the source and noticed that 
no ice absorption features are present, while the 10 fim region indicate the combination of sili- 
cate emission and absorption. Since our ISOPHOT-S and Spitzer/IRS spectra are consistent 
with this flnding, we decided not to assign any type to this source. All the ISOCAM-CVF, 
ISOPHOT-S, and Spitzer spectra indicate signiflcant flux changes and changes in the spectral 
slope. 



WXCha is a TTauri-type binary with a separation of 0'.'79 and P.A. of 55° (iGhez et al. 



19971 ). Both the ISO and Spitzer spectra contain all the flux from both components, thus, 
they can be safely compared. 



CV Cha is a single T Tauri-type star (iGuenther et al.ll2007l : iMelol 120031 1 . although it shows 
variable radial velocity ( Reipurth et al. 20021 ). The Spitzer/IRS slit contained only CVCha, 
while the ISOPH OT-S beam included a lso CWCha (at a separation of 11"). Resolved IRAC 
photometry from iLuhman et al.l (120081 ) indicates that the flux ratio of CW Cha to CV Cha 
is 0.12-0.16 between 3.6 and 8^m. This, and the fact that CWCha falls at the very edge of 
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the ISOPHOT-S beam ensures that CW Cha has neghgible contribution to the ISOPHOT- 
S spectrum. Thus, we conclude that both the ISOPHOT-S and the Spitzer/IRS spectra 
represent the flux from CV Cha only, thus they can be compared. 



HD 98800 (also known as TV Crt) is a quadruple T Tauri system. It consists of a binary 
with a projected separation of C'S (or 38 AU), P. A. 0°, and each component is a spectroscopic 
binary with a separation of about 1 AU (jPrato et al.ll2001t iBoden et al.ll2005l ). Both the ISO 
aperture and the IRS slit c ontain all four obj ects, so they can be compared. The IR excess is 
attributed to HD 98800 B JPrato et al.lboOli l. HD 988 00 B displays no ex cess emission below 
about 5.5 /im, indicating an inner cleared-out region (IFurlan et al.l 120071 ). 



HD 100453 is an isolated binary Herbig star. The primary has a spectral type of A9 Ve. 
The companion, HD 100453 B, is situated at a s eparation of y'OQ ( or 120 AU), and is 110 
times fainter in the Kg band than the primary (IHabart et al.ll2006[). The compan ion is a 
M4.0-M4.5V star with no significant excess in the K or L band (I Collins et al.l I20091). t hus it 
has no contribution to the plotted ISOPHOT-S or Spitzer/IRS spectra. iMeeus et al.l (|2002[ ) 
explain the lack of the 10 /im silicate feature with the depletion of small, hot silicate grains in 
the disk of HD 100453 (a factor of 100-500 depletion compared to AB Aur is needed in their 
model to be able to consiste nt with the observed IS O-SWS spectrum). Using higher signal- 
to-noise ratio TIMMI2 data. Ivan Boekel et al.l (|2005[ ) confirmed the lack of the 10 /im silicate 
emission, and attributed it to grain growth. The comparison of our ISOPHOT-S and IRS 
spectra reveal some differences in the 8-12 /xm wavelength range, while the flux at A < 8 //m 
is constant. It is possible that this variability is due to a weak, variable silicate emission 
feature (Fig. IT^ . Confirmation of this idea would require a detailed modeling of the IRS 
spectrum, which has the highest signal-to-noise ratio among the 10 //m spectra published so 
far for HD 100453. The star shows no strong accretion activity and is photometrically stable 
( ICoUins et al.l 120091 : iGarcia Lopez et al.ll2006l : iGuimaraes et al.ll2006l ). 



HD 10423 7 is a spectroscop ic binary consisting of a Herbig A4V primary and a K3-type 
companion (IBohm et al.l 120041 ). Based on X-ray observations, iTesta et al.l ( 120081 ) identified 
four more lo w-mass stars in the v icinity, HD 104237 B-E, with separations between 1'.'365 and 
14'.'88, while iGrady et al.l (l2004j ) found five low-mass sources in their optical, near-IR and 
mid-IR images. In addition to the Herbig Ae star itself, HD 104237-6 (corresponding to the 
X-ray source HD 104237 E) and HD 104237-2 have IR excesses, but the H erbig star HD 104237 
dominates the mid-IR integrated light of the region ( IGrady et al.ll2004l ). While HD 104237- 
6 falls outside of both the ISOPHOT-S beam and the Spitzer/IRS slit, HD 104237-2 is 
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included in both. The Spitzer sht is centered between the Herbig star and HD 104237-2. 
Offset correction was apphed using the position of the Herbig star. Although the observed 
Spitzer/IRS spectrum may contain 75-100% of the flux of HD 104237-2, its contribution 
to the much brighter primary is negligible, thus we conclude that the ISOPHOT-S and 
Spitzer/IRS spectra can be compared. Nevertheless, the ISOPHOT-S spectra alone already 
shows significant variability (see also Fig. [T6|l . 



DKCha is a deeply embedde d active Herbig Ae star. It displays brightness variations of 
more than Imag in the K band (IHughes et al.lll99l[ ). The source is redder when fainter, and 
the amplitude of the brightness fluctuations decreases with increasing wavelength. They 
explain the observed fl ux changes w i th a combination of variable circumstellar extinction 
and intrinsic variability. iPrusti et al.l ( 119941 ) found no evidence for extended mid-IR emission 
(with mult i- diaphragm observations between 3.8 and 18.6 /xm, meaning that the emitting area 
is less than 5.4"). We note th at there is a significant discrep ancy between the ISOPHOT-S 
(this work) and the ISO-SWS (lAcke &: van den Anckerl 120041 ) spectra of the source, pointing 
to possible variability. The Spitzer data points are close to saturation and are not very 
reliable. 



HD 135344 B is an F4 Ve-type Herbig star fivan Boekel et al.ll2005h . lOoucet et all (120061 ) 

resolved the source at 20.5 /zm, and modeled the emission with a disk with an outer radius 
of 200 AU ( o r 1"4 which is within both the ISOPHOT-S and the Spitzer/IRS aperture). 
Grady et al.l (120091 ) observed significant changes between 1-10 /im, both in fiux and in the 
shape of the SED. At wavelengths shortward of 4.75 /xm, at least 13% variations are present, 
while at 10 /um, they found 60% variability. 



HD 139614 is an Herbig Ae star (IDunkin et al.lll997h . lYudin et al.l (119991 ) detected po- 
larimetric variability on a time-scale of days, which they explain in terms of a model of 
a pole-on disk with dust clouds or comet-like bodies moving in a circumstellar envelope 
around the star. We found no evidence in the literature that the source might be extended 
at mid-IR wavelengths, thus we conclude that the ISOPHOT-S and the Spitzer spectra can 
be safely compared. The two spectra agree within the measurement uncertainties, more- 
over, the TIMMI2 spectrum of Ivan Boekel et al.l (120051 ) is also identical with them. Thus, 
HD 139614 is c onstant at mid-IR w avelengths (nor does it show any significant variability in 
the optical, see lMeeus et al.lll998l ). 
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HD 141569 is a B9.5/A0 Herbig star (iJaschek fc Jascheklll992l:lDunkin et al.lll997h. Coro- 



nagraphic o bservations of t he sc attered light by lAugereau et al.l ( Il999l ). [Weinberger et aL 



( 119991 ) ■ and iMouillet et al.l (12001) revealed that its disk is asymmetric, and consists of two 
rings with a gap between them. [Weinberger et al.l (120001 ) claims t hat the disk is compose d 
of secondary debris material, although it still contains some gas (jZuckerman et al.lll995l ). 
HD 141569 has two M-type companions: HD 141569 B (M2V) and HD 141569 C (M4V) at 
project ed distances of 7^.^6 ( A-B) and I'H (B-C); K -band flux ratios are: 5.3 (A/B) and 11.9 
(A/C) dPirzkal et al]ll997l : Iweinberger et al.lboooh . The Spitzer/IRS slit included only the 
primary (A), while the ISOPHOT-S aperture included all three components. Due to their 
faintness and distance from component A, components B and C have negligible contribution 
to the ISOPHOT-S flux, thus the ISOPHOT-S and Spitzer/IRS fluxes can be compared. 



HD 142666 is a Herbig Ae star (Ivan Boekel et alll2005h. Mid-IR spect ra of HD 142666 



obtain ed w ith different instruments are published in ISylvester et al.l (119971 ) ; iBouwman et al. 



(l200ll) , and Ivan Boekel et al.l ( 120051 ) . HD 142666 displa ys large photometric variatioris (iMalfait et al, 
19981 give a visual amplitude of about 1.2 mag, while Lecavelier Des Etangs et al. 20051 de- 
tected 0.3 mag variability with Stromgren photometry). The star appears redder when 
fainter, which they interpret as variable dust extinction in a close to edge-on disk with 
clumpy dust distribution. Since beam effects do not play a role here, we can safely compare 
the ISOPHOT-S and Spitzer spectra. 



HD 142527 is an isolated Fe-type Herbig star (Ivan Boekel et al.l l2005l . and references 



therein). Based on spectro-astrometric observations. iBaines et al 



( 120061 ) consider it a possi- 



ble subarcsecond binary. The circums tellar environment of the star was resol ved at various 
wave lengths (see e.g. near-IR i mages in IPukagawa et al" mid-IR images in lFujiwara et al 



20061 . or millimeter images in lVerhoeff et al.l 1201 ll ). Based on these observations, the star is 
surrounded by an inner disk, a n outer disk, an d a ha lo arou nd the inner disk regio ns. Spa- 
tially resolved 18 /um images in iFujiwara et al.l ( l2006l ) and in IVerhoeff et al.l (120111 ) indicate 
that the emission is coming from within about 1" of the star, which is fully included in both 
the ISOPHOT-S beam and the Spitzer/IRS slit. Thus, the spectra can be compared, and 



they i ndicate significant variability. An additional 8-13 /jm spectrum from Ivan Boekel et al. 
( 120051 ) also confirms this (Fig. [T7|) . 



EX Lup is the prototype of EXors, a class of young stars sh owing sporadic optical out- 
burst separated by extended quiescent periods (iHerbia 119771 ). The outbursts are pow - 
ered by enhanced accretion from the circumstellar disk onto the star (IHerbig et al.l 120011 ). 
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EXLup itself brightened by AKsiSmag in 1955-56 and again in 2008, and produced sev- 
eral smaller, 1-3 mag brightenings in-between. EXLup was extensively studied during its 
outburst in 2008, a nd the observation s indicate that it brightened in the whole 0.4-100 yum 
wavelen gth range (IJuhasz et al.l |2012[ ) . By comparing pre-outburst and outburst Spitzer 
spectra, Abraham et al. ( 20091 ) observed a significant change in the shape of the 10 /jm sil- 
icate feature: while before the outburst, the feature indicated mostly amorphous silicates, 
during outburst, se veral narrow pe a ks an d shoulders appeared, indicating the presence of 
crystalline silicates. [Abraham et al.l ( l2009l ) explain the appearance of the crystalline silicates 
by annealing in the surface layer of the inner disk by heat from the outburst. In Fig. [2l we 
plot ISOPHOT-S spectra from 1997 and Spitzer spectra from 2004 and 2005. These are all 
quiescent spectra, when the visual magnitude of EXLup was between 12.5 mag and 14 mag. 



HD 144432 is a pre-main sequence binary consisting of HD 144432 A, an A9/F0-type Her- 
big star primary and HD 144432 B, a KO-5-type T Tauri secondary, with a separation of 1.4", 
P.A. 0° dPerez et al.ll2004i Istelzer et allbood : IPukagawa et~aDl201oh. The secondary is a fac- 
tor o f 9-14 fainter than the primary in the K and H bands (jPerez et al.ll2004j : iFukagawa et al. 
20101). The secondary displays phot ospheric emission until K band, and the IR excess is at- 
tributed exclusively to the primary (jPerez et al.ll2004j ). Both the ISOPHOT-S aperture and 
the Spitzer/IRS slit include both components and the observed spectra represent the total 
fiux of the two objects (the S pitzer slit is oriented along the binary, north-south). Keck 
segment-tilting observations by iMonnier et al.l (l2009l ) resolved the source at 10.7 //m, which 
is Ri40mas. There is a ~10% difference between the ISOPHOT-S and the Spitzer/IRS 
spectra plotted in Fig. [21 Severa l othe r mid-IR spectra are available on HD 144432 in the 
literature, see e.g. [Sylvester et al.l ( Il996[ ) ; [Leinert et al.l (|2004| ) ; Ivan Boekel et al.l (120051 ) . The 
first of these is consistent with the Spitzer/IRS fiux level, the second one is consistent with 
the ISOPHOT-S fiux level, and the third one is lower than the others by another ^^10%, 
supporting our finding that the source is variable in the mid-IR regime. 



HR5999 is a bin ary star consis ting of a Herbig Ae primary, and a T Tauri secondary 
called Rossiter 3930 (lRossiteiill955l). The separation is 1"4, the P.A. is 111°, and the K -band 



fiux ratio is between 20-35 ( IStecklum et al.lll995 



Ghez et al. 



1997 



pri mary itself may be a sp e ctrosc opic binary ( Tjin A Djie et al. 19891 ). although observations 



Leinert et al.lll997h . The 



by ICorporon fc Lagrangd (119991 ) do not confirm this . HR5999 is a kn own optical variable, 
whose "bursts" are attributed to unsteady accretion (jPerez et al.lll992l) . Multi-epoch BVRI- 
JHKLMNQ band photometry by lHutchinson et al.l (Il994l ) indicated that the object is indeed 
variable at o ptical wavelengths but not mu ch above 4.8 /im. An ISOCAM CVF spectrum 
published in lAcke fc van den Anckerl ( 120041 ) is consistent with our ISOPHOT-S spectrum 
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until 8.5 /iin, but is lower by about 20% at l onger wavelengths, indicating a possibly variable 
silicate feature. ISiebenmorgen et al. J2000h noted that HR5999 is a point source at these 
wavelengths. Using VLTI/MIDI, iPreibisch et al.l (120061 ) resolved the mid-IR emission from 
HR5999 and they give sizes between 5 and 25 mas (1-5 AU), depending on the fitted model. 



WL 16 is a single B8-A7 Herbig star (jSimon et al. I I1995I : iRatzka et ahlboOslJLuhman fc Riekel 
I999I ). Images at different mid-IR wavelength s between 7.9 and 24.5 /xm (both in the con- 
tinuum and in the PAH bands) published by iRessler fc Barsonyl (120031 ) revealed that the 



central star is surrounded by a bright disk with a size of 7"x3'!5 



diameter of 900 AU), viewed at an inclination of 62.2°. iGeers et al. 



corresponding to a disk 
(I2OO7I ) also resolve both 



the 8-13 fim continuum and the 8.6, 11.2, and 12.7 urn PAH fea tures, but give smaller radial 
extents (50-60 AU, or 0'.'4-0'.'5) than iRessler &: Barsonyl (120031 ). In their study, the 3.3 /im 
feature is undetected, indicating the dominance of ionized PAHs. Our ISOPHOT-S spec- 
trum of WL 16 in Fig. 1 agrees well with the ISO-S WS spectrum, a nd th e slight differences 
between these and the 10 //m spectrum published in iHanner et al.l (Il992l ) are probably due 
to the different beam sizes of the instruments. 



MWC 863 is a pre-main sequence binary, consist ing of a Herbig Ae pri mary, HD 150193 A, 
and an F9-type TTauri secondary, HD 150193 B (ICarmona et al.l 120071 ). Their separation 



is I'.'l, the P.A. is 227° (iReipurth fc Zinnecker 
K and L' bands, respectiv ely Jzinnecker et al ' 



ponent has no I R excess (iFukagawa et al.l 



graphic images, iFukagawa et al 



I993I), flux ratios a re 6-8 and 4-17 in the 



1991 



200 



Koreskol |2002| ) . The secondary com- 

^, and references therein). Using corona- 



(120031 ) resolved the disk around HD 150193 A in the H 
(190 AU). During Keck segment-tilting observations by 
Monnier et al.l (120091 ) the source remained unresolved at 10.7 yum, which gives an upper 



band, and gave a ra dius of 1'.'3 



limit of FW HM=35 mas for the emitt i ng re gion. We note that while the ISO-SWS spectrum 
pubhshed in Acke fc van den Anckeil (2004) is identical with ou t ISOPHOT-S spectrum, the 
UKIRT/CGS3 spectrum taken by lSylvester fc Mannings! (I2OOOI ) is significantly different both 
in shape an in flux level, pointing to possible changes in the composition and size distribution 
of the silicate dust grains. 



AK Sco is a double-line spectroscopic binary consisting of approximately equal mass F5- 
type stars with a proj ected separation of 0.14 AU (~1 mas); the two stars are surrounded by 
a circumbinary disk (lAlencar et al.l l2003l ) . Both the ISO and Spitzer spectra include both 
components. Its optical light curve shows a "roughly constant maxi mum light level beiri g 
interrupted at irregular intervals by minima about a magnitude deep" (I Andersen et al .11 19891 ) . 
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Our ISOPHOT-S and S pitzer/IRS spectra agree within the measurement uncertainties, and 
a TIMMI2 spectrum by iPrzygodda et al.l (120031 ) also has the same flux level. 



MWC865 (also known as y921Sco and CD-4 2°11721) is a BOIV-type star embedded 
in a small, thick dark cloud (IBoersma et al.l l2009l and references therein). It illuminates a 
reflection nebula ( van den Bergh fc HerbstI 19751 ). I ts distance is quite uncertain (values be - 
tween 400 and 2500 pc appear in the literature, e.g. iBoersma et al.ll2009l: iLopes et al.lll992l). 



and there is also some debate about its pre-main sequence natu re (iBorges Fernandes et al. 



20071 ). although it is generally considered to be a H erbig Be star. 



Ha 



jart et al.l (120031 ) found 



four IR companions within 13" of MWC 865, while IWang fc LoonevI (l2007t ) found about 10 
YSO candidate s in th e vicinity (within about 70") of MWC 865 using Spitzer/IRAC images. 
Boersma et al.l ( 120091 ) shows Spitzer/IRAC images between 3.6 and 8.0 fim and TIMMI2 
images at 11.9 fim of MWC 865 and its surroundings. These images reveal two bright peaks, 
one coinciding with the star, the other one being a close-by patch 4"3 to the north. There is 
an arc stretching towards the south-east as well, at a distance of 11" from the star. TIMMI2 
long-slit spectra indicate that the central star ha s a featureless spectru m, while the close-by 
patch and the arc display strong PAH emission. IBoersma et al.l (120091 ) combine the spectra 
of these components and compare it with an ISO-SWS spectrum. They find that the ISO- 
SWS spectrum shows a similar shape, but has higher flux levels due to the larger aperture 
of ISO-SWS and the extent of the PAH emission. Our ISOPHOT-S spectrum of MWC 865 
also has a very similar shape to the ISO-SWS spectrum, but even higher flux levels. We 
conclude that these differences are also due to the difference in the apertures and their exact 
position with respect to the PAH-emitting patches. We note that significant offset correction 
was applied to the ISOPHOT-S spectrum of MWC 865, introducing extra uncertainty in the 
absolute flux levels. 



51 Qph is a single Herbig star ( iBaines et al.l 120061 ) . which is photometrically constant 



(iLecavelier Pes Etangs et al 



20051). The source appeared point-like in the 18 /im images 



of IJayawardhana et al.l (1200 ll) . givin g an upper limit of 0"5 (65 AU) for the radius of the 
emitting region. iLeinert et al.l (120041 ) resolved the source at 12.5 /im with VLTI/MIDI and 
gave a half-light radius of 7mas (0.5 AU). Our conclusion is that the ISOPHOT-S and the 
Spitzer/IRS spectra can be compared and the observed yariabi hty is real. We n ote that 



ground-based mid-IR spectra published in ISylvester et al.l (Il996[ ) and iLeinert et al.l (120041 ) 
are consistent with the Spitzer/IRS flux level. 



-42 - 



HD 163296 is a well-known single ( ICorporon &: Lagrangelll999l: iPirzkal et al.l 119971) Her- 
big Ae star. Using optical and near-IR coronagrapliic images, Fukagawa et al. (j2010 ) and 



Grady et al.l (1200(1) coul d detect the disk in scattered light out to a radius of 3'.'6 (440 AU). 



Wisniewski et al.l (l2008l ) found that the optical scattered light is va riable in time due to 
variable self-shadowing. During Keck segment-tilting observations by iMonnier et al.l (120091 ) 
the source remained un resolved at 10 . 7 /xm, which gives an upper limit of FWHM=35mas 
for the emitting region. iDoucet et al.l ( l2006l ) resolved the disk thermal emissi on at 20.5 iim 



and d erived an outer disk radius of 1'.'6 (200 AU). The source was resolved by lLeinert et al. 



(120041 ) using VLTl/MlDl data, and the half-light r adius at 12.5 urn is 7mas (0 .8 AU). The 
source appeared point-like in the 18 /im images of IJayawardhana et al.l (l200ll), givi ng an 
upper limit of O'.'S (60 AU) for the radius of the emitting region. ISitko et al.l (120081 ) con- 
ducted a long-term 3—13 /xm spectroscopic monitoring of HD 163296. Most of their data 
are consistent with no variability larger than 10%, except at one epoch when the 1— 5/im 
flux significantly increased, accompanied by a slight increase in the flux of the 10 /im silicate 
feature. Although no Spitzer/IRS low resolution spectra are available for this source, our 
ISOPHOT-S spectrum is consistent with an ISO-SWS spectr um (taken on the same day), 
and also consistent with the TIMMI2 spectrum published i n Ivan Boekel et al.l (l2005r) . ob- 
tained in March 2003, and the MIDI spectrum published in iLeinert et al.l (l2004l), obtained 
in June 2003, while it is 10-40% higher than that obtained by iKessler-Silacci et al. ( 2005 ) 
with the Keck/LWS in August 1999 and in June 2000. 



HD 169142 is a Herbig Ae star (ISylvester et al.lll996l ). High spatial resolution mid-IR 
observations at 3.3, 10.8, and 1 8.2 nm show that the emission comes from a d isk with an outer 
radius of 150 AU (1") at most (iHabart et al.ll2006l : iJavawardhana et al.ll2001l). SED modelin g 



gave disk radii in the range of 100-300 AU (0'.'7-2") (iDominik et al 



20031: 



Dent et al.ll2006r ). 



HD 169 142 has a possible companion, 2MASS J18242929-2946559, located at a distance 
of 9'.'3 (iGrady et al.ll2007l ). The 2MASS source is a binary weak-line TTauri with 130 mas 
separation, and is 11 times fainter than HD 169142 in the K band. The Spitzer/IRS slit 
contains only HD 169142, while the ISOPHOT-S beam includes both components. However, 
the contribution of the 2MASS source at mid-IR wavelengths is probably negligible, thus the 
Spitzer and the ISO spectra can be compared. The spectra plotted in Fig. 1 agree within the 
measurement uncertainties, thus we conclude that this sour ce is constant in the 5 -11.5 
regime. The flux levels o f the UKIRT/CGS3 spe ctrum from ISylvester et al.l (Il996[ ) and the 
TIMMI2 spectrum from Ivan Boekel et al.l (120051 ) also agree with our spectra until ~9 fim, 
while they differ by up to 50% above 9 iim (F ig. ITSl) . The source is photometrically stable 
at optical wavelengths (iGuimaraes et al.ll2006l ). 
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VV Ser is a single (IPirzkal et al. 



1997; 



type optical variability ( Herbst fc Shevchenkol 



Leinert et a. 



1999 



19971) Herbig star. It d isplays UXor- 



Rostopchina et al.ll200l[ ). i.e. 1-3 mag 



deep minima in th e V band due to var iable obscuration by dust clumps in a nearly edge-on 
circumstellar disk. Habart et al. found no other source emitting at 10 ^m in the 1' x 1' 



vicinity of V VSer. 



Li et al. 



(Il99. 



) found no extended emission around the source in JHK 



filters, while iLeinert et al.l ( 1200 lb) found faint extended emission with a FWHM of 0.6" in 



JHK with speckle interferometry. lEisner et al.l ( l2004j ) resolved the source at 2.2 /im with the 
PTI, an d determined a size of 1.5 -4.5 mas (depending on the model used) for the emitting 
region. iPontoppidan et al.l ( 120071 ) discovered nebulous mid-IR emission extending over 4' 



centered on VVSer. The nebulosity is due to transiently heated dust grains. Although we 
cannot exclude the possibility that this extended emission has different contributions to the 
ISOPHOT-S and the Spitzer/IRS spectra, the fact that the Spitzer/IRS flux level is higher 
points to real variability. 



OO Ser is a deeply emb edded YSO that prod uced an outburst in 1995 and became 4.6 mag 
brighter in the K ban d (IHodapp et al.l Il996l ). reaching its maximum brightness in 1995 
October (lHodapplll999l). The eruption caused brightening at a wide wavelength range from 
2.2 /xm to 100 ;um (I Kospal et al.l 120071 ). After peak brightness the object started a gradual 
fading. The ISOPHOT-S spectra plotted in Fig. |2] cover a period of 20 months, and the first 
one was obtained 4 months after the max i mum brightness occurred. The gradual fading, 
which is discussed in detail in lKospal et al.l ( 120071 ) . is well visible in the figure. The source is 
surrounded by a small elongated nebula, approxir aately 15" in diameter, best seen in the K 
band, but also discernible at 3.6, 4.5 and 5.8 //m (IHodapp et al.lll996l : iKospal et al.ll2007D . 



CK 1 is an embedded protostellar binary with a separation of 1'.'5, and P.A. of 10° (lEiroa fc Leinert 



19871 ). The southern component is brighter than the northern one, and t heir flux ratio de- 
creases from 4.2 at 0. 9 jim to 3.3 at 3.5 /x m to 1.1-2.8 between 8-12.9 jjm (lEiroa et al.lll987l : 
Eiroa fc Leinertlll987l : ICiardi et al.ll2005l ). Both the ISO and the Spitzer spectra represent 
the sum of the spectra of CK 1 north and CK 1 south (the Spitzer s l it is o riented along the 
separation of the binary), thus, they can be compared. ICiardi et al.l (120051 ) present resolved 
8-13 ^m spectra of the two components, and the sum of the two spectra as well. This sum 
agrees very well with our ISOPHOT-S and Spitzer/IRS data. 



[SVS7 6] Ser 4 is a small, dens e cluster of deeply embedded low- and int ermediate mass 
YSOs (jPontoppidan et aL 2004). Judging fro m the JHK photome try of Eiroa fc Casali 
(119891 ) and the Spitzer 3.6-24.0 //m fluxes from iHarvey et al.l (120071 ). most of the flux in 
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the Spitzer slit comes from SVS4/9 and SVS4/10, while in addition to these two sources, 
SVS4/4, SVS4/8, SVS4/11 and SVS4/5 also contribute to the ISO flux. SVS4/5 is especially 
bright above 5 fim (very red source). Source confusion explains the differences between the 
ISO and Spitzer spectra. We note that in the case of the ISOPHOT-S spectrum, predicted 
background was subtracted, making the absolute level of the measured flux probably more 
uncertain than 10%. 



CK 2 is an embedded YSO in the Serpens star-forming region. We found no indication 
of extended mid-IR emission in the literature. The source is close to the sensitivity limit of 
ISOPHOT-S, resulting in a very low signal-to-noise, especially above 6 /im. The ISOPHOT-S 
spectrum also suffers from memory effect, and a predicted background is subtracted, which, 
considering how faint the source is, introduces additional uncertainties. Th us, we do not an- 



alyze the variability of CK2. Nevertheless, the ISOCAM CVF spectrum of [Alexander et al. 



( 120031 ) is significantly lower than our Spitzer/IRS spectrum, indicating possible variability. 



S CrA is a PMS binary with a separa tion of r.'3-r.'5, P . A of 147°-160° K-band flux ratio 



1995; 


Prato et al. 


2003 


)• 



[Baier et al. 1985: Zinnecker et al. 1991: Chelli et al. 



era 



for the two components and determined a spectral type of G5 Ve for the primary (S CrA NW) 
and K5Ve for the secondary (S Cr A SE). S CrA is a well-known o ptical variable, with V-band 
variations of up to ^2 mag (e.g. iKardopolov fc Filipevlll98ll ). (jGrahamlll992l ) interpreted 
the photometric and spectroscopic changes of the object as variable accretion onto the star. 
Both the ISOPHOT-S aperture and the Spitzer/IRS slit included the whole flux of both 
components, thus, they can be compared. 



R CrA, T CrA, and HH 100 IRS are part of Coronet, a small cluster of YSOs. The two 
brightest sources, R CrA and T CrA, are Herbig stars, while t he remaining dozen ( HH 100 IRS 
among others) are deeply embedded low-mass protostars ( iForbrich et al.l 120071 . and refer- 
ences therein). RCrA illuminates the reflection n ebula NGC6729. Both R CrA and the 
nebula are highly variable at optical wavelengths ( Graham fc Phillips 1987h. HH 100 IRS 



is also highly variable, with amplitudes up to 2.5 mag in the K band ( lAxon et al.l Il982 



Reipurth & Wamste 


cer 


1983) 


naries (iTakami et al. 


20031). I 



19831 ) . Both R CrA and T CrA a re possible spectro-ast rometric bi- 



mid-IR emission larger than ^5" around R CrA and T CrA. Interestingly R CrA displays the 
4.27 ixm feature of solid CO2 ice in absorption and the 10 ^m silicate feature in emission at 
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the same time ( iNummelin et al.ll200ll : lAcke &: van den Anckerll2004[ ). No Spitzer spectra are 
available for RCrA and TCrA, but t he comparison of our ISOPHO T-S spectra in Fig. [2] 
and the ISO-SWS spectra published in lAcke &: van den Anckerl (120041 ) indicate very similar 
flux levels. For HH 100 IRS, variability can be studied by comparing the two ISOPHOT-S 
spectra. 



VV CrA is a TTauri-type binary with a separation of 2'.'10, P.A. of 44°, and K-band flux 
ratio of 8.6. The optical se condary becomes b righter around the J band, and dominates the 
flux at longer wavelengths (IChelli et al. II 19 951 ). iRatzka et al.l (120081 ) report that the secondary 
faded significantly and "it is now fainter than the primary even in the N-band". They also 
show 10 fim MIDI spectra for each component: VVCrANE shows deep silicate absorption, 
while VVCrASW is featureless. Our ISO spectrum contains flux from both components, 
while the Spitzer/IRS spectrum in centered on the NE component, but contains significant 
contribution from the SW as well. For this reason, we do not analyze variability for this 
source. 



WWVul is an isolated Herbig Ae star, which shows UXor-like fadings at optic al wave- 
lengths, although its disk probably has a larger inclination than for other UXors (jPugach 
20041 : iKozlova et aPbooeh . WWVul was observed twice with ISOPHOT-S and once with 
Spitzer. Both ISOPHOT-S spectra are offset-corrected, but the large correction apphed for 
the earlier one (TDT: 17600465) makes this more uncertain than the 10% assigned as a final 
uncertainty of the ISOPHOT-S spectra. For this reason, this spectrum in not used in our 
variability analysis, and we base our conclusions on the comparison of the other ISOPHOT-S 
spectrum (TDT: 51300108) and the Spitzer measurement. 



BP - |-40 4124 and LkHa224 are members of a cluster of young stars (iHenning et al. 



19981 . and references therein). The brightest members of the cluster are the B2-type BD -|-40 4124 
and the B5-type LkHa224, but several fainter FGKM stars can be found within a few ar- 
cminutes; some of these are included in the ISOPHOT-S aperture. No Spitzer spectrum is 
obtained for these sources, however an ISO-SWS is available. Comparison of our ISOPHOT-S 
spectrum with the ISO-SWS reveals no significant differences. BD + 40 4124 displays small- 
scale (<0.4 mag) quasi -periodic variability at optical wavelengths (IHerbst fc Shevchenko 
1999l : lstrom et al.lll972h . 



PV Cep is an embedded YSO surrounded by the cometary reflection ne bula RNO 125 
(iThe et al.l Il994j ). Both the star and the nebula are strongly variable (e.g. IScarrott et al. 
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199ll : iLorenzetti et al.ll201l[ ). iKun et al.l (1201 if ) analyzed the brightness changes of PVCep 



in the whole optical-infrared regime, and found significant flux variations below 10 fim. They 
conclude that the reason for the variability is a combination of changing accretion rate and 
changing circumstellar extinction. We found no indication in the literature that the object 
is extended at mid-IR wavelengths. The fact that the Spitzer spectrum has a higher flux 
than the ISOPHOT-S points to real variability. 



SV Cep is spectro-astrometric binary Herbig star ( IWheelwright et al.l l2010l ). Its opti 
cal light curve is characteristic of UXor-type variables with irregular, q uasi-periodic varia- 



tions d ue to inhomogeneitie s in a nearly edge-on circumstellar disk (e.g. iRostopchina et al. 
2000). lJuhasz et al.l ( 120071 ) investigated the long-term 0.55-100 /um variability of SVCep 
and found a correlation between the optical and far-IR light curves. In the mid-IR regime, 
they found lower variability amplitudes (a peak-to-peak variability of ~40% at 10 fim) and a 
hint for anti-correlation with the optical changes. They interpret the variability using a self- 
shadowed disk and a puffed- up inner rim with variable inner rim height. Their model also 
contained an optically thin envelope to account for the constant 25 //m flux. Our ISOPHOT- 
S spectrum plotted in Fig. |2] and the ISO-SWS spectrum obtained 9 months previously 
( Acke &: van den Ancker 20041 ) are consistent with each other. 



REFERENCES 

Abraham, P., Kospal, A., Csizmadia, S., Kun, M., Moor, A., & Prusti, T. 2004, A&A, 428, 
89 

Abraham, P., et al. 2009, Nature, 459, 224 

Acke, B., & van den Ancker, M. E. 2004, A&A, 426, 151 

Acosta-Pulido, J. A., & Abraham, P. 2003, in ESA Special Publication, Vol. 481, The 
Calibration Legacy of the ISO Mission, ed. L. Metcalfe, A. Salama, S. B. Peschke, & 
M. F. Kessler, 95 

Alencar, S. H. P., Johns-Krull, C. M., & Basri, G. 2001, AJ, 122, 3335 

Alencar, S. H. P., Melo, C. H. F., Dullemond, C. P., Andersen, J., Batalha, C, Vaz, L. P. R., 
& Mathieu, R. D. 2003, A&A, 409, 1037 

Alexander, R. D., Casali, M. M., Andre, P., Persi, P., & Eiroa, C. 2003, A&A, 401, 613 

Andersen, J., Lindgren, H., Hazen, M. L., & Mayor, M. 1989, A&A, 219, 142 



-47- 

Arellano Ferro, A., & Giridhar, S. 2003, A&A, 408, L29 

Augereau, J. C, Lagrange, A. M., Mouillet, D., & Menard, F. 1999, A&A, 350, L51 

Axon, D. J., Allen, D. A., Bailey, J., Hough, J. H., Ward, M. J., & Jameson, R. F. 1982, 
MNRAS, 200, 239 

Baier, G., Keller, E., Weigelt, G., Bastian, U., & Mundt, R. 1985, A&A, 153, 278 
Bailey, J. 1998, MNRAS, 301, 161 

Baines, D., Oudmaijer, R. D., Porter, J. M., & Pozzo, M. 2006, MNRAS, 367, 737 

Barsony, M., Ressler, M. E., & Marsh, K. A. 2005, ApJ, 630, 381 

Bary, J. S., Leisenring, J. M., & Skrutskie, M. F. 2009, ApJ, 706, L168 

Beck, T. L., Schaefer, G. H., Simon, M., Prate, L., Stoesz, J. A., & Howell, R. R. 2004, ApJ, 
614, 235 

Beichman, C. A., et al. 1984, ApJ, 278, L45 
Boden, A. F., et al. 2005, ApJ, 635, 442 

Boersma, C., Peeters, E., Martin- Hernandez, N. L., van der Wolk, G., Verhoeff, A. P., 
Tielens, A. G. G. M., Waters, L. B. F. M., & Pel, J. W. 2009, A&A, 502, 175 

Bohm, T., Catala, C., Balona, L., & Carter, B. 2004, A&A, 427, 907 

Bontemps, S., et al. 2001, A&A, 372, 173 

Borges Fernandes, M., Kraus, M., Lorenz Martins, S., & de Araiijo, F. X. 2007, MNRAS, 
377, 1343 

Bouwman, J., Meeus, G., de Koter, A., Hony, S., Dominik, C., & Waters, L. B. F. M. 2001, 
A&A, 375, 950 

Bowey, J. E., & Adamson, A. J. 2001, MNRAS, 320, 131 

Brandeker, A., Liseau, R., Artymowicz, P., & Jayawardhana, R. 2001, ApJ, 561, L199 
Bregman, J., Rank, D., Sandford, S. A., & Temi, P. 1993, ApJ, 410, 668 
Brinch, C., Crapsi, A., Hogerheijde, M. R., & J0rgensen, J. K. 2007, A&A, 461, 1037 
Carmona, A., van den Ancker, M. E., & Henning, T. 2007, A&A, 464, 687 



-48 - 

Casali, M. M. 1991, MNRAS, 248, 229 

Chelli, A., Cruz-Gonzalez, I., & Reipurth, B. 1995, A&AS, 114, 135 
Chiang, E. I., & Goldreich, R 1997, ApJ, 490, 368 

Ciardi, D. R., Telesco, C. M., Packham, C, Gomez Martin, C, Radomski, J. T., De Buizer, 
J. M., Phillips, C. J., & Barker, D. E. 2005, ApJ, 629, 897 

Cohen, M. 2003, in ESA Special Publication, Vol. 481, The Calibration Legacy of the ISO 
Mission, ed. L. Metcalfe, A. Salama, S. B. Peschke, & M. F. Kessler, 135 

Collins, K. A., et al. 2009, ApJ, 697, 557 

Connelley, M. S., Reipurth, B., & Tokunaga, A. T. 2008, AJ, 135, 2496 
Corporon, P., & Lagrange, A. 1999, A&AS, 136, 429 

Correia, S., Zinnecker, H., Ratzka, T., & Sterzik, M. F. 2006, A&A, 459, 909 

Cutri, R. M., et al. 2003, VizieR Onhne Data Catalog, 2246, 

Davies, J. K., Evans, A., Bode, M. P., & Whittet, D. C. B. 1990, MNRAS, 247, 517 

Decin, L., Morris, P. W., Appleton, P. N., Charmandaris, V., Armus, L., & Houck, J. R. 
2004, ApJS, 154, 408 

Dent, W. R. P., Torrelles, J. M., Osorio, M., Calvet, N., & Anglada, G. 2006, MNRAS, 365, 
1283 

Desert, P., Boulanger, P., & Puget, J. L. 1990, A&A, 237, 215 
Dewarf, L. E., & Dyck, H. M. 1993, AJ, 105, 2211 

Dominik, C, DuUemond, C. P., Waters, L. B. P. M., & Walch, S. 2003, A&A, 398, 607 

Doucet, C, Habart, E., Pantin, E., Dullemond, C, Lagage, P. O., Pinte, C, Duchene, G., 
& Menard, F. 2007, A&A, 470, 625 

Doucet, C, Pantin, E., Lagage, P. O., & Dullemond, C. P. 2006, A&A, 460, 117 

Duchene, G., Bcust, H., Adjali, P., Konopacky, Q. M., & Chez, A. M. 2006, A&A, 457, L9 

Dullemond, C. P., Dominik, C, & Natta, A. 2001, ApJ, 560, 957 



-49 - 



DuUemond, C. P., Henning, T., Visser, R., Geers, V. C, van Dishoeck, E. F., & Pontoppidan, 
K. M. 2007, A&A, 473, 457 

Dunkin, S. K., Barlow, M. J., & Ryan, S. G. 1997, MNRAS, 290, 165 

Eaton, N. L., & Herbst, W. 1994, in Bulletin of the American Astronomical Society, Vol. 26, 
American Astronomical Society Meeting Abstracts #184, 933 

Eiroa, C., & Casali, M. M. 1989, A&A, 223, L17 

Eiroa, C., & Leinert, C. 1987, A&A, 188, 46 

Eiroa, C., Lenzen, R., Leinert, C., & Hodapp, K. 1987, A&A, 179, 171 

Eisner, J. A., Lane, B. P., Hillenbrand, L. A., Akeson, R. L., & Sargent, A. 1. 2004, ApJ, 
613, 1049 

Espaillat, C, Furlan, E., D'Alessio, P., Sargent, B., Nagel, E., Calvet, N., Watson, D. M., 
& MuzeroUe, J. 2011, ApJ, 728, 49 

Feigelson, E. D., & Kriss, G. A. 1989, ApJ, 338, 262 

Flaherty, K. M., MuzeroUe, J., Rieke, G., Gutermuth, R., Balog, Z., Herbst, W., Megeath, 
S. T., & Kun, M. 2011, ApJ, 732, 83 

Forbrich, J., et al. 2007, A&A, 464, 1003 

Fujiwara, H., et al. 2006, ApJ, 644, L133 

Fukagawa, M., Tamura, M., Itoh, Y., Hayashi, S. S., & Oasa, Y. 2003, ApJ, 590, L49 

Fukagawa, M., Tamura, M., Itoh, Y., Kudo, T., Imaeda, Y., Oasa, Y., Hayashi, S. S., & 
Hayashi, M. 2006, ApJ, 636, L153 

Fukagawa, M., et al. 2010, PASJ, 62, 347 

Furlan, E., et al. 2006, ApJS, 165, 568 

— . 2007, ApJ, 664, 1176 

Gabriel, C., Acosta-Pulido, J., Heinrichsen, 1., Morris, H.. & Tai, W. 1997, in Astronom- 
ical Society of the Pacific Conference Series, Vol. 125, Astronomical Data Analysis 
Software and Systems VI, ed. G. Hunt & H. Payne, 108-111 

Garcia Lopez, R., Natta, A., Testi, L., & Habart, E. 2006, A&A, 459, 837 



Geers, V. C, van Dishoeck, E. F., Visser, R., Pontoppidan, K. M., Augereau, J., Habart, 
E., & Lagrange, A. M. 2007, A&A, 476, 279 

Geers, V. C, et al. 2006, A&A, 459, 545 

Ghez, A. M., McCarthy, D. W., Patience, J. L., & Beck, T. L. 1997, ApJ, 481, 378 

Gotz, W. 1980, Information Bulletin on Variable Stars, 1747, 1 

Grady, C. A., et al. 2000, ApJ, 544, 895 

— . 2004, ApJ, 608, 809 

— . 2007, ApJ, 665, 1391 

— . 2009, ApJ, 699, 1822 

Graham, J. A. 1992, PASP, 104, 479 

Graham, J. A., & Phillips, A. C. 1987, PASP, 99, 91 

Grinin, V. P., Shakhovskoi, D. N., Shenavrin, V. I., Rostopchina, A. N., & Tambovtseva, 
L. V. 2002, Astronomy Reports, 46, 646 

Guenther, E. W., Esposito, M., Mundt, R., Covino, E., Alcala, J. M., Cusano, F., & Steck- 
lum, B. 2007, A&A, 467, 1147 

Guimaraes, M. M., Alencar, S. H. P., Corradi, W. J. B., & Vieira, S. L. A. 2006, A&A, 457, 
581 

Giirtler, J., Schreyer, K., Henning, T., Lemke, D., & Pfau, W. 1999, A&A, 346, 205 

Habart, E., Natta, A., & Kriigel, E. 2004, A&A, 427, 179 

Habart, E., Natta, A., Testi, L., & Carbillet, M. 2006, A&A, 449, 1067 

Habart, E., Testi, L., Natta, A., & Vanzi, L. 2003, A&A, 400, 575 

Haisch, Jr., K. E., Barsony, M., Ressler, M. E., & Greene, T. P. 2006, A J, 132, 2675 

Haisch, Jr., K. E., Greene, T. P., Barsony, M., & Stabler, S. W. 2004, AJ, 127, 1747 

Hammersley, P. L., & Jourdain de Muizon, M. 2003, in ESA Special Publication, Vol. 481, 
The Calibration Legacy of the ISO Mission, ed. L. Metcalfe, A. Salama, S. B. Peschke, 
& M. F. Kessler, 129 



-51 - 

Banner, M. S., Tokunaga, A. T., & Geballe, T. R. 1992, ApJ, 395, Llll 

Harvey, P., Merm, B., Huard, T. L., RebuU, L. M., Chapman, N., Evans, II, N. J., & Myers, 
R C. 2007, ApJ, 663, 1149 

Henning, T. 2010, ARA&A, 48, 21 

Henning, T., Burkert, A., Launhardt, R., Leinert, C, & Stecklum, B. 1998, A&A, 336, 565 
Herbig, G. H. 1977, ApJ, 217, 693 
— . 1990, ApJ, 360, 639 

Herbig, G. H., Aspin, C., Gilmore, A. C., ImhofF, C. L., & Jones, A. F. 2001, PASP, 113, 
1547 

Herbst, W., Herbst, D. K., Grossman, E. J., & Weinstein, D. 1994, AJ, 108, 1906 
Herbst, W., & Shevchenko, V. S. 1999, AJ, 118, 1043 
Hessman, F. V., & Guenther, E. W. 1997, A&A, 321, 497 
Hirose, S., & Turner, N. J. 2011, ApJ, 732, L30 

Hodapp, K., Hora, J. L., Rayner, J. T., Pickles, A. J., & Ladd, E. F. 1996, ApJ, 468, 861 
Hodapp, K. W. 1999, AJ, 118, 1338 

Houck, J. R., et al. 2004, in Presented at the Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference, Vol. 5487, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, ed. J. C. Mather, 62-76 

Hughes, J. D., Hartigan, P., Graham, J. A., Emerson, J. P., & Marang, F. 1991, AJ, 101, 
1013 

Hutchinson, M. G., Albinson, J. S., Barrett, P., Davies, J. K., Evans, A., Goldsmith, M. J., 
& Maddison, R. C. 1994, A&A, 285, 883 

Isobe, S., Norimoto, Y., & Kitamura, T. 1988, PAS J, 40, 89 

Jaschek, C, & Jaschek, M. 1992, A&AS, 95, 535 

Jayawardhana, R., Fisher, R. S., Telesco, C. M., Pina, R. K., Barrado y Navascues, D., 
Hartmann, L. W., & Fazio, G. G. 2001, AJ, 122, 2047 



-52 - 

Jensen, E. L. N., Dhital, S., Stassun, K. G., Patience, J., Herbst, W., Walter, F. M., Simon, 
M., & Basri, G. 2007, AJ, 134, 241 

Joblin, C., Tielens, A. G. G. M., Geballe, T. R., & Wooden, D. H. 1996, ApJ, 460, L119 

Juhasz, A., Prusti, T., Abraham, P., & DuUemond, C. P. 2007, MNRAS, 374, 1242 

Juhasz, A., et al. 2010, ApJ, 721, 431 

— . 2012, ApJ, 744, 118 

Kardopolov, V. I., & Filipev, G. K. 1981, Soviet Ast., 25, 457 
Keller, L. D., et al. 2008, ApJ, 684, 411 
Kenyon, S. J., et al. 1994, AJ, 107, 2153 

Kessler-Silacci, J. E., Hillenbrand, L. A., Blake, G. A., & Meyer, M. R. 2005, ApJ, 622, 404 
Kohler, R., Ratzka, T., Herbst, T. M., & Kasper, M. 2008, A&A, 482, 929 
Koresko, C. D. 2002, AJ, 124, 1082 

Kospal, A., Abraham, P., Prusti, T., Acosta-Pulido, J., Hony, S., Moor, A., & Siebenmorgen, 
R. 2007, A&A, 470, 211 

Kospal, A., Salter, D. M., Hogerheijde, M. R., Moor, A., & Blake, G. A. 2011, A&A, 527, 
A96 

Kospal, A., et al. 2008, MNRAS, 383, 1015 
Kospal, A.., et al. 2012, A&A 

Kozlova, O. v., Shakhovskoi, D. N., Rostopchina, A. N., & Alekseev, I. Y. 2006, Astro- 
physics, 49, 151 

Kun, M., et al. 2011, MNRAS, 413, 2689 

Lafreniere, D., Jayawardhana, R., Brandeker, A., Ahmic, M., & van Kerkwijk, M. H. 2008, 
ApJ, 683, 844 

Lagage, P., et al. 2006, Science, 314, 621 



-53- 



Laureijs, R. J., Klaas, U., Richards, R J., Schulz, B., & Abraham, P. 2003, The ISO Hand- 
book, Volume IV - PHT - The Imaging Photo-Polarimeter, ed. Kessler, M. P., Mueller, 
T. G., Leech, K., Arviset, C, Garcia-Lario, P., Metcalfe, L., Pollock, A., Prusti, T., 
& Salama, A. 

Lecavelier Des Etangs, A., Nitschelm, C., Olsen, E. H., Vidal-Madjar, A., & Ferlet, R. 2005, 
A&A, 439, 571 

Leinert, C., Abraham, P., Acosta-Pulido, J., Lemke, D., & Siebenmorgen, R. 2002, A&A, 
393, 1073 

Leinert, C., Beck, T. L., Ligori, S., Simon, M., Woitas, J., & Howell, R. R. 2001a, A&A, 
369, 215 

Leinert, C., & Haas, M. 1989, ApJ, 342, L39 

Leinert, C., Haas, M., Abraham, P., & Richichi, A. 2001b, A&A, 375, 927 

Leinert, C., Haas, M., Mundt, R., Richichi, A., & Zinnecker, H. 1991, A&A, 250, 407 

Leinert, C., Richichi, A., & Haas, M. 1997, A&A, 318, 472 

Leinert, C., Zinnecker, H., Weitzel, N., Christou, J., Ridgway, S. T., Jameson, R., Haas, M., 
& Lenzen, R. 1993, A&A, 278, 129 

Leinert, C., et al. 2004, A&A, 423, 537 

Li, W., Evans, II, N. J., Harvey, P. M., & Colome, C. 1994, ApJ, 433, 199 
Lim, J., & Takakuwa, S. 2006, ApJ, 653, 425 
Liu, M. C., et al. 1996, ApJ, 461, 334 

Lopes, D. P., Damineh Neto, A., & de Preitas Pacheco, J. A. 1992, A&A, 261, 482 

Lorenzetti, D., Larionov, V. M., Giannini, T., Arkharov, A. A., Antoniucci, S., Nisini, B., 
& Di Paola, A. 2009, ApJ, 693, 1056 

Lorenzetti, D., et al. 2011, ApJ, 732, 69 

Luhman, K. L., Allen, P. R., Espaillat, C., Hartmann, L., & Calvet, N. 2010, ApJS, 186, 
111 

Luhman, K. L., & Rieke, G. H. 1999, ApJ, 525, 440 



-54- 

Luhman, K. L., et al. 2008, ApJ, 675, 1375 

Magakian, T. Y., & Movsesian, T. A. 2001, Astrophysics, 44, 419 
Malfait, K., Bogaert, E., & Waelkens, C. 1998, A&A, 331, 211 
Manfroid, J., ct al. 1991, A&AS, 87, 481 

McCabc, C, Ghez, A. M., Prato, L., Duchene, G., Fisher, R. S., & Telesco, G. 2006, ApJ, 
636, 932 

Meeus, G., Bouwman, J., Dominik, G., Waters, L. B. F. M., & de Koter, A. 2002, A&A, 
392, 1039 

Meeus, G., Waelkens, G., & Malfait, K. 1998, A&A, 329, 131 

Meeus, G., Waters, L. B. F. M., Bouwman, J., van den Ancker, M. E., Waelkens, G., & 
Malfait, K. 2001, A&A, 365, 476 

Melo, G. H. F. 2003, A&A, 410, 269 

Monnier, J. D., Tuthill, P. G., Ireland, M., Gohen, R., Tannirkulam, A., & Perrin, M. D. 
2009, ApJ, 700, 491 

Moore, T. J. T., & Emerson, J. P. 1992, MNRAS, 259, 381 

— . 1994, MNRAS, 271, 243 

Morales-Galderon, M., et al. 2009, ApJ, 702, 1507 
— . 2011, ApJ, 733, 50 

Mouillet, D., Lagrange, A. M., Augereau, J. G., & Menard, F. 2001, A&A, 372, L61 

Murakawa, K., Gya, S., Pyo, T., & Ishii, M. 2008, A&A, 492, 731 

MuzeroUe, J., Galvet, N., Hartmann, L., & D'Alessio, P. 2003, ApJ, 597, L149 

MuzeroUe, J., Luhman, K. L., Briceno, G., Hartmann, L., & Galvet, N. 2005a, ApJ, 625, 906 

MuzeroUe, J., Megeath, S. T., Flaherty, K. M., Gordon, K. D., Rieke, G. H., Young, E. T., 
& Lada, G. J. 2005b, ApJ, 620, L107 

MuzeroUe, J., et al. 2009, ApJ, 704, L15 



-55 - 

Myers, P. C, Puller, G. A., Mathieu, R. D., Beichman, C. A., Benson, P. J., Schild, R. E., 
& Emerson, J. P. 1987, ApJ, 319, 340 

Nagel, E., D'Alessio, P., Calvet, N., Espaillat, C, Sargent, B., Hernandez, J., & Forrest, 
W. J. 2010, ApJ, 708, 38 

Natta, A. 1993, ApJ, 412, 761 

Nummelin, A., Whittet, D. C. B., Gibb, E. L., Gerakines, P. A., & Chiar, J. E. 2001, ApJ, 
558, 185 

Osorio, M., D'Alessio, P., MuzeroUe, J., Calvet, N., & Hartmann, L. 2003, ApJ, 586, 1148 
Perez, M. R., van den Ancker, M. E., de Winter, D., & Bopp, B. W. 2004, A&A, 416, 647 
Perez, M. R., Webb, J. R., & The, P. S. 1992, A&A, 257, 209 
Pirzkal, N., Spillar, E. J., & Dyck, H. M. 1997, ApJ, 481, 392 
Pojmanski, G. 1997, Acta Astron., 47, 467 

Pontoppidan, K. M., Dullemond, C. P., Blake, G. A., Evans, II, N. J., Geers, V. C., Harvey, 
P. M., & Spiesman, W. 2007, ApJ, 656, 991 

Pontoppidan, K. M., van Dishoeck, E. P., & Dartois, E. 2004, A&A, 426, 925 

Prato, L., Greene, T. P., & Simon, M. 2003, ApJ, 584, 853 

Prato, L., Simon, M., Mazeh, T., Zucker, S., & McLean, I. S. 2002, ApJ, 579, L99 
Prato, L., et al. 2001, ApJ, 549, 590 

Preibisch, T., Kraus, S., Driebe, T., van Boekel, R., & Weigelt, G. 2006, A&A, 458, 235 

Prusti, T., Clark, F. O., Whittet, D. C. B., Laureijs, R. J., & Zhang, C. Y. 1991, MNRAS, 
251, 303 

Prusti, T., Natta, A., & Palla, F. 1994, A&A, 292, 593 

Przygodda, F., van Boekel, R., Abraham, P., Melnikov, S. Y., Waters, L. B. F. M., & Leinert, 
C. 2003, A&A, 412, L43 

Pugach, A. F. 2004, Astronomy Reports, 48, 470 

Ratzka, T., Kohler, R., & Leinert, C. 2005, A&A, 437, 611 



-56- 



Ratzka, T., Leinert, C, Przygodda, F., & Wolf, S. 2008, in The Power of Optical/IR Interfer- 
ometry: Recent Scientific Results and 2nd Generation, ed. A. Richichi, F. Delplancke, 
F. Paresce, & A. Chelli, 519-522 

Ratzka, T., et al. 2009, A&A, 502, 623 

RebuU, L. M. 2011, ArXiv e-prints 

Reipurth, B. 1985, A&AS, 61, 319 

Reipurth, B., & Bally, J. 1986, Nature, 320, 336 

Reipurth, B., Lindgren, H., Mayor, M., Mermilliod, J., & Cramer, N. 2002, AJ, 124, 2813 
Reipurth, B., & Wamsteker, W. 1983, A&A, 119, 14 
Reipurth, B., & Zinnecker, H. 1993, A&A, 278, 81 
Ressler, M. E., & Barsony, M. 2003, ApJ, 584, 832 

Richichi, A., Leinert, C, Jameson, R., & Zinnecker, H. 1994, A&A, 287, 145 
Roche, P. F., Aitken, D. K., & Smith, C. H. 1994, MNRAS, 269, 649 
Rossiter, R. A. 1955, Publications of Michigan Observatory, 11, 1 

Rostopchina, A. N., Grinin, V. P., & Shakhovskoi, D. N. 2001, Astronomy Reports, 45, 51 

Rostopchina, A. N., Grinin, V. P., Shakhovskoi, D. N., The, P. S., & Minikulov, N. K. 2000, 
Astronomy Reports, 44, 365 

Sandell, G., & Weintraub, D. A. 2001, ApJS, 134, 115 

Sargent, B. A., et al. 2009, ApJS, 182, 477 

Scarrott, S. M., Rolph, C. D., & Tadhunter, C. N. 1991, MNRAS, 249, 131 
Scarrott, S. M., & Wolstencroft, R. D. 1988, MNRAS, 231, 1019 

Schmidt, T. O. B., Neuhauser, R., Seifahrt, A., Vogt, N., Bedalov, A., Helling, C., Witte, 
S., & Hauschildt, P. H. 2008, A&A, 491, 311 

Shiba, H., Sato, S., Yamashita, T., Kobayashi, Y., & Takami, H. 1993, ApJS, 89, 299 

Siebenmorgen, R., Natta, A., Kruegel, E., & Prusti, T. 1998, A&A, 339, 134 



-57- 

Siebenmorgen, R., Prusti, T., Natta, A., & Miiller, T. G. 2000, A&A, 361, 258 
Simon, M., et al. 1995, ApJ, 443, 625 

Sipos, N., Kospal, A., Grossman, E. J., & Weinstein, D. 2011, A&A, in prep 
Sitko, M. L., et al. 2008, ApJ, 678, 1070 

Skemer, A. J., Close, L. M., Hinz, P. M., Hoffmann, W. F., Kenworthy, M. A., & Miller, 
D. L. 2008, ApJ, 676, 1082 

Sloan, G. C., Hayward, T. L., AUamandola, L. J., Bregman, J. D., Devito, B., & Hudgins, 
D. M. 1999, ApJ, 513, L65 

Smith, K. W., Lewis, G. F., Bonnell, I. A., Bunclark, P. S., & Emerson, J. P. 1999, MNRAS, 
304, 367 

Stapelfeldt, K. R., et al. 1995, ApJ, 449, 888 

Stecklum, B., Eckart, A., Henning, T., & Loewe, M. 1995, A&A, 296, 463 

Stelzer, B., Robrade, J., Schmitt, J. H. M. M., & Bouvier, J. 2009, A&A, 493, 1109 

Sterzik, M. F., Melo, C. H. F., Tokovinin, A. A., & van der Bliek, N. 2005, A&A, 434, 671 

Straizys, V., Corbally, C. J., Kazlauskas, A., & Cernis, K. 2002, Baltic Astronomy, 11, 261 

Strom, K. M., Strom, S. E., Breger, M., Brooke, A. L., Yost, J., Grasdalen, G., & Carrasco, 
L. 1972, ApJ, 173, L65 

Strom, S. E., Vrba, F. J., & Strom, K. M. 1976, AJ, 81, 314 

Sylvester, R. J., & Mannings, V. 2000, MNRAS, 313, 73 

Sylvester, R. J., Skinner, C. J., & Barlow, M. J. 1997, MNRAS, 289, 831 

Sylvester, R. J., Skinner, C. J., Barlow, M. J., & Mannings, V. 1996, MNRAS, 279, 915 

Takami, M., Bailey, J., & Chrysostomou, A. 2003, A&A, 397, 675 

Testa, P., Huenemoerder, D. P., Scliulz, N. S., & Ishibashi, K. 2008, ApJ, 687, 579 

The, P. S., de Winter, D., & Perez, M. R. 1994, A&AS, 104, 315 

Tielens, A. G. G. M. 2008, ARA&A, 46, 289 



-58- 



Tjin A Djie, H. R. E., The, P. S., Andersen, J., Nordstrom, B., Finkenzeller, U., & Jankovics, 
I. 1989, A&AS, 78, 1 

Torres, G. 2004, AJ, 127, 1187 

Turner, N. J., Carballido, A., & Sano, T. 2010, ApJ, 708, 188 

van Boekel, R., Juhasz, A., Henning, T., Kohler, R., Ratzka, T., Herbst, T., Bouwman, J., 
& Kley, W. 2010, A&A, 517, A16 

van Bockcl, R., Min, M., Waters, L. B. F. M., de Koter, A., Dominik, C, van den Ancker, 
M. E., & Bouwman, J. 2005, A&A, 437, 189 

van Boekcl, R., Waters, L. B. F. M., Dominik, C, DuUemond, C. R, Tielens, A. G. G. M., 
& de Koter, A. 2004, A&A, 418, 177 

van den Bergh, S., & Herbst, W. 1975, AJ, 80, 208 

van Dishoeck, E. F. 2004, ARA&A, 42, 119 

Verhoeff, A. P., et al. 2011, A&A, 528, A91 

Vinkovic, D., Ivezic, Z., Jurkic, T., & Elitzur, M. 2006, ApJ, 636, 348 
Vinkovic, D., & Jurkic, T. 2007, ApJ, 658, 462 
Wang, S., & Looney, L. W. 2007, ApJ, 659, 1360 
Watson, D. M., et al. 2009, ApJS, 180, 84 

Weinberger, A. J., Becklin, E. E., Schneider, G., Smith, B. A., Lowrance, P. J., Silverstone, 
M. D., Zuckerman, B., & Terrile, R. J. 1999, ApJ, 525, L53 

Weinberger, A. J., Rich, R. M., Beckhn, E. E., Zuckerman, B., & Matthews, K. 2000, ApJ, 
544, 937 

Wheelwright, H. E., Oudmaijer, R. D., & Goodwin, S. P. 2010, MNRAS, 401, 1199 

Wisniewski, J. P., Clampin, M., Grady, C. A., Ardila, D. R., Ford, H. C., Golimowski, D. A., 
niingworth, G. D., & Krist, J. E. 2008, ApJ, 682, 548 

Wooden, D. H., Bell, K. R., Barker, D. E., & Woodward, C. E. 2000, in Bulletin of the 
American Astronomical Society, Vol. 32, American Astronomical Society Meeting 
Abstracts, 1482 



-59- 



Yudin, R. V., Clarke, D., & Smith, R. A. 1999, A&A, 345, 547 

Zinnecker, H., Christou, J. C, Ridgway, S. T., & Probst, R. G. 1991, in Astronomical 
Society of the Pacific Conference Series, Vol. 14, Astronomical Society of the Pacific 
Conference Series, ed. R. Elston, 270-272 

Zuckerman, B., Forveille, T., & Kastner, J. H. 1995, Nature, 373, 494 



This preprint was prepared with the A AS macros v5.2. 



Table 1. Catalog of spectral observations of young stellar objects obtained with ISO and Spitzer 
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Table 1 — Continued 



Coordinates ISO Spitzer 
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08-Aug-2004 


em 




em 


sil. em. 


yes 
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em 
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18 
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47800913 


B 
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em 
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34300318 




24-Oct-1996 
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ll-Mar-1996 
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FO 


19 
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BO 


06- Apr- 1996 
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L 


K7 
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Table 1 — Continued 



Coordinates ISO Spitzcr 



Name 


Mass 


Sp.typc 


RA(2000) 


Doc(2000) 


TDT 


Corr. 


Obs. date 


AOR 


Corr. 


Obs. date 


Si-O 


Ice 


PAH 


Typo 


Var. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


WW Vul 


I 


A3 


19 25 58.74 


+21 12 31.4 


17600465 


MBO! 


ll-May-1996 


16828672 




20-.Iun-2006 


cm 






sil. cm. 


yes 












51300108 


BO 


12-Apr-1997 


















BD +40 4124 


I 


B3 


20 20 28.25 


+41 21 51.5 


15900568 


B 


24- Apr- 1996 












cm 


PAH dom. 




LkHa 224 


I 


F9 


20 20 29.36 


+41 21 28.5 


14201271 


BO 


07-Apr-1996 












cm 


PAH dom. 




PV Cep 


L 


GO-KO" 


20 45 53.94 


+67 57 38.7 


14302274 


BO 


08-Apr-1996 


12287488 




23-Oct-2004 


abs 


abs 




sil. abs. 


yes 


SV Cep 


I 


AO 


22 21 33.20 


+73 40 27.1 


56201203 


B 


31-May-1997 








cm 






sil. cm. 





Note. — For a detailed description of the columns of this table, see Sec. [3] 
^''■^ISO observations for LDN 1489 IRS and OO Ser were taken in mapping mode. 
^^■^ Coordinates for Ccd 111 IRS 5 come from a Spitzer/IRAC 3.6 /./.m image. 

^^^2MASS coordinates for the member of the SVS4 cluster that is brightest at >2.2 ^im, SVS4/9 ([EC92] 95). 
^"^^ Coordinates for S Cr A is from lPrato et al.| i2003l "). 

SIMBAD gives a spectral type of A5, but the star is probably a low-mass object of GO-KO iMagak ian &: Movsesia ^ |2001H . 
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Fig. 1. — Spectra of the spectrophotometric standard star HD 181597 (HR7341). In the up- 
per panel, sohd hnes show the observations (TDT 53102202 for ISOPHOT-S, AOR 28709632 
for Spitzer) , while dotted and dashed hnes show the template and model spectra. The mid- 
dle panel displays the ratio between the observed and the Cohen template spectrum, while 
the lower panel is the ratio between the observed and the Decin model spectrum. 
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Fig. 2. — ISOPHOT-S and Spitzer/IRS spectra for our object. Where clearly detected, 
vertical dashed lines mark the wavelengths of molecular ice bands and vertical dotted lines 
indicate the PAH features. 
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Fig. 2. — (Continued) 
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Fig. 2. — (Continued) 
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Fig. 2. — (Continued) 
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Fig. 2. — (Continued) 
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Fig. 2. — (Continued) 
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Fig. 3. — Ratios of spectra obtained at different epochs. Where several spectra existed, we 
took the two most extremes. Error bars represent the quadratic sum of the errors of each 
spectrum. 
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Fig. 3. — (Continued) 
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Fig. 3. — (Continued) 
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Fig. 4. — Histograms showing the distribution of magnitude changes for our sources, (a)- 
(d): Magnitude changes in the continuum; (e)-(h): Magnitude changes in the feature (PAH 
features for PAH dominated objects, sihcate feature for objects showing sihcate emission 
or absorption). Different hues are plotted for magnitude changes calculated between the 
ISOPHOT-S and Spitzer/IRS spectra (long-term variabihty) and for magnitude changes 
calculated for spectra taken with the same instrument (short-term variability) . 
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Fig. 5. — Synthetic magnitude changes. The graphs display the brightness changes of our 
targets in the continuum and in the PAH or silicate band. The typical uncertainty of 0.1 mag 
is represented by the dark gray circle in the middle of each panel. The vertical light gray 
stripe marks the area where sources with constant continuum but variable silicate feature 
fall. The diagonal light gray stripe marks the location of sources whose continuum and 
feature changed with similar amplitudes. 
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Fig. 6. — Example for a cosmic glitch which was only partly corrected by the standard PIA 
deglitching algorithm, producing a spike in the spectrum (HD 97300, TDT 62501316, Pixel 
60). In our interactive processing all data points after t=330s are manually discarded. 
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Fig. 7. — Example for memory effect during the observation of CRCha (TDT 6250f703). 
The measured signal of Pixel 93 (A=8.45 ^m) decreases rather than increases with time, due 
to recent illumination of the pixel by a bright source. The measured data points lie above 
the expected transient curve (solid line), leading to an overestimation of the derived flux. 




Fig. 8.— Dark current observation before the observation of CRCha (TDT 62501703). The 
observed signal (thick sohd hne) is well above the mission average value (dashed line, with 
la errors indicated by dotted lines). 
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Fig. 9. — Contribution of the zodiacal background (dashed hne) to the measured flux of 
BF Ori (solid line). Since no dedicated background measurement was performed, the zodiacal 
spectrum was estimated from the DIRBE 4.9 and 12 /xm data points extracted for the position 
and date of the BF Ori observation. 
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Fig. 10. — Example for spectral artifacts introduced by off-center location of a point source 
in the ISOPHOT-S beam. The object (Cedll21RS4, TDT 62501217) was offset by +4.8" 
and —0.6" in satellite Y- and Z-direction, respectively. Vertical marks show the wavelengths 
where the most obvious artifacts, in the forms of discontinuities, appear. 
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Fig. 11. — Upper panel: Spectrum of the standard star HD 144284 (dashed hne, TDT 
52401702) compared with a stellar template prediction (solid line). Lower panel: Ratio 
of the measured flux values to the stellar template. 
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Fig. 12. — Ratio of measured flux values to stellar template prediction for 3 wavelengths 
(all normal star observations are included). Uncertainty values representative of the whole 
ensemble are computed separately at high flux level (expressed in %) and at low flux level 
(in mJy). 
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Fig. 13. — Typical photometric uncertainties for the ISOPHOT-S wavelength range. Mul- 
tiplicative errors (upper panel) are valid for brighter sources, while additive errors (lower 
panel) are for fainter sources (<1 Jy). 
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Fig. 14.— Mid-IR spectra of HD 100453. The ISOPHOT-S (black line) and Spitzer/IRS 
(light gray line) s pectr a are from this work, the T1MM12 (light gray) spectrum is from 
van Boekel et al.l ( 120051 ). The bottom panel shows the difference between the spectra ob- 
tained in 2003 and 2004. 
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Fig. 15.— Mid-IR spectra of HD 169142. The ISOPHOT-S (dashed hne) and Spitzer/IRS 
(black hne) spectra ar e from this work, the TIMMI2 (hght gray) spe ctrum is from 
van Boekel et al.l ((20051), the UKIRT/CGS3 (dark gray) spectrum is from ISylvester et al. 
Jl996h . 
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Fig. 16.— Mid-IR spectra of HD 104237. The ISOPHOT-S (dotted and dashed lines) and 
Spitzer/IRS (black line ) spectra are from this work, the TIMMI2 (gray) spectrum is from 
van Boekel et al.l koO^ . 
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Fig. 17.— Mid-IR spectra of HD 142527. The ISOPHOT-S (dashed hne) and Spitzer/IRS 
(black line) spectra are from this work, the T1MM12 (gray) spectrum is from Ivan Boekel et al. 



